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Abstract 

Dramatic progress has occurred in the last two decades in understanding the physical 
processes and events leading up to, and transpiring during the eruption of a classical nova. 
The mechanism whereby a white dwarf accreting hydrogen-rich matter from a low-mass 
main-sequence companion produces a nova eruption has been understood since 1970. The 
mass-transferring binary stellar configuration leads inexorably to thermonuclear runaways 
on the white dwarf surfaces, mass ejection, and photon emission so luminous it can be 
detected at distances of megaparsecs. This review summarizes the efforts of many researchers 
in understanding the physical processes which generate nova eruptions; the effects upon 
nova eruptions of different binary-system parameters (e.g. chemical composition or mass of 
the white dwarf, different mass accretion rates); the possible metamorphosis from dwarf to 
classical novae and back again; and observational diagnostics of novae, including X-ray and 
7 -ray emission, and the characteristics and distributions of novae in globular clusters and 
in extragalactic systems. While the thermonuclear-runaway model remains the successful 
cornerstone of nova simulations, it is now clear that a wide variety of physical processes, and 
three-dimensional hydrodynamic simulations, will be needed to explain the rich spectrum of 
behavior observed in erupting novae. 
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I. INTRODUCTION AND SCOPE OF THE REVIEW 

The first systematic observations of erupting novae were made by Chinese imperial as- 
trologers two thousand years ago. Their underlying reason fcr observing novae was the same 
as that of modern astronomers. Anything that brightens overnight by fifteen magnitudes is 
interesting; however the use to which sky watchers put novae has changed in the last two 
millenia. Instead of being regarded as astrological omens, novae are now acknowledged to 
be powerful cosmological tools. Novae are excellent standard candles; they are extremely 
luminous at maximum, and their rates of decline are tightly correlated with their absolute 
magnitudes at maximum (Arp, 1956). The better a standard candle is understood, the less 
chance there is for misusing or misinterpreting it. This alone justifies the detailed study of 
classical novae. For many workers in the field, however, the rich potpourri of theoretical 
physics, numerical simulation, and range of observational phenomena involving novae keeps 
research going at a brisk pace regardless of these systems’ cosmological value. 

It is ironic that extragalactic novae were at first part of a serious stumbling block to 
cosmology. The supernova of 1885 in the nucleus of M31 was incorrectly argued to be a 
classical nova by Shapley in his famous debate with Curtis. Only with the resolution of 
Cepheids in M31 by Hubble were the relative brightnesses and distinct natures of classical 
novae and supernovae finally understood by astronomers. Had the supernova of 1885 not 
thrown astronomers off the track, novae might have helped settle the nature of the nebulae 
years earlier. 

A wealth of observational data on classical novae was summarized by Payne-Gaposchkin 
(1957) and by McLaughlin (1960). It was clear to these authors that classical novae ejected 
significant amounts of matter at velocities of hundreds or thousands of kilometers per second 
into space. Faint, hot blue stars remained after nova eruptions, and the stars appeared about 
as bright after eruption as before. The intrinsically most luminous novae were observed to 
belong to the fastest “speed class” *.e., to fade most rapidly. Very fast nova decline by ~3 
magnitudes in ~10 days; very slow novae may require a year to dim by the same amount. 

The crucial breakthrough in understanding why novae erupt was the observation by 
Walker (1954) and by Kraft (1959, 1964) that novae occur in white dwarf-main sequence 
binary systems where hydrogen-rich matter is being accreted onto the white dwarf. Ac- 
cretion occurs via a disk which is generated by the red dwarf overflowing its Roche Lobe. 
Considerable energy and angular momentum losses allow the disk to feed matter onto the 
white dwarf’s surface. 

Starrfield, Truran, Sparks and Kutter (1972) and Starrfield, Sparks and Truran (1974a, 
19746) demonstrated convincingly that hydrogen-rich envelopes on the surfaces of white 
dwarfs undergo thermonuclear runaways (TNRs). These explosions are violent enough to 
eject matter at high velocity if the envelopes are enriched in CNO nuclei relative to the Sun. 
Bath and Shaviv (1976) showed that radiation pressure can drive winds for extensive periods 
from novae. Subsequent observations in the infrared, ultraviolet, and X-ray wave bands 
have demonstrated that the bolometric luminosities of novae remain remarkably constant 
at the Eddington luminosity for months after the visual luminosity decreases. The physical 
processes leading to this observed behaviour include large bolometric corrections created by 
the thinning of the wind (Bath and Shaviv 1976), the contraction and heating of the remnant 
envelope (Prialnik, Shara and Shaviv 1978) and dust formation and heating (Bode 1983). 

This appealingly straightforward scenario of nova eruptions was summarized a decade 
ago in an extensive review by Gallagher and Starrfield (1978). Since then, many researchers 
have been exploring the wide range of physical processes involved in nova eruptions, and the 
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dramatic effects these have on the subsequent eruption. The aim of this review is to describe 
the efforts of these workers, over the last ten years, to progress from a simple understanding 
of the basic nova mechanism to a detailed picture of the three-dimensional, time-dependent, 
highly nonlinear processes that govern nova eruptions. 

Several valuable reviews have been published in recent years on subjects related to cat- 
aclysmic binaries. An extraordinary compilation of finder charts, accurate coordinates, and 
outburst properties of Galactic novae has been assembled and published by Duerbeck (1987). 
This remarkable atlas of ~200 classical novae is of enormous value to all researchers in the 
field. A further, important compilation of cataclysmic binary data and interpretation is 
that of Patterson (1984). Specialized reviews have been published on the subjects of cat- 
aclysmic binaries with moderate magnetic fields attached to the white dwarf (the DQ Her 
stars: Warner 1985); the strongly magnetic AM Herculis stars (Liebert and Stockman 1985); 
the evolution of magnetic cataclysmics (Lamb and Melia 1987); IUE ultraviolet studies of 
galactic novae (Starrfield 1986); dust in nova systems (Bode and Evans 1985); and nucle- 
osynthesis in novae and the effects of /3 + -unstable nuclei and their effects on nova eruptions 
(Truran 1985). I will touch only lightly on these topics. 

The ~2-5 magnitude outbursts of dwarf novae, repeating on time scales of weeks to years, 
are seen only in low mass transfer-rate binaries. They probably occur when a recurring disk 
instability rapidly causes much of the accretion disk to be dumped onto the white dwarf. 
The significant release of gravitational potential energy manifests itself as an erupting dwarf 
nova. Dwarf novae are briefly considered as metamorphic forms of classical nova systems 
within the context of the hibernation theory (§VI.l). Finally, the evolution of binary stars 
to and from the cataclysmic stage, and the reason for the total lack of cataclysmics in the 
2-3 hour period range will be discussed in a forthcoming review by Livio (1989). 


II. THE PRE-OUTBURST STAGE 
III RADIAL ACCRETION OF HYDROGEN-RICH MATERIAL 
ONTO A WHITE DWARF 

The first generation of TNR models (Starrfield, Sparks and Truran 1974a; 1974b; Pri- 
alnik, Shara and Shaviv 1978, 1979) were based upon hydrogen-rich envelopes already in 
place on the white dwarf. While these models succeeded in reproducing the energetics and 
early light curves observed in classical novae, it was clear that the assumed envelope masses 
and thermal structures were not necessarily self consistent. Numerical simulations of the 
envelope accretion and evolution were first carried out by Kutter and Sparks (1980), Nariai, 
Nomoto, and Sugimoto (1980), and by Prialnik et al. (1982). 

The most important parameter in determing whether a TNR occurs was found to be 
the mass accretion rate M. Rapid mass accretion (M> 10 - ®M©yr -1 ) produces significant 
compressional heating at the hydrogen-rich envelope base. A mild TNR occurs because 
the accreted mass is not very degenerate, and the envelope puffs up into a red giant-like 
configuration. Only for M ^ 10~®,M©i/r -1 is there sufficient degeneracy built up at the 
envelope base prior to ignition to result in a strong TNR. The recent inclusion of boundary 
layer heating in accretion calculations (Shaviv and Starrfield 1988; Prialnik, Kovetz and 
Shara 1988) implies even a lower value, M & 10” 10 M©yr -1 for a strong TNR (see §11.4). 
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At first sight, these upper limits conflict with observations of “old” novae (systems 10-100 
years after eruption). Patterson (1984) has compiled observations of old-nova brightnesses 
and distances. He finds that M~ 10 is implied by the observed brightnesses of 
most old novae. Yet if the white dwarfs in old novae continue accreting matter from their 
companions at 10~ 8 Mgyr -1 , simulations indicate that they will not again erupt at classical 
novae. This is contrary to the findings of Ford (1978) who showed conclusively that classical 
novae must suffer hundreds or thousands of TNRs during their lifetimes This is because ~50 
novae/year are observed to erupt in M31 (see §VII.2). Over a Hubble time (~ 1.5 x 10 10 
years) roughly 7.5 x 10 11 nova eruptions occur in M31. This is comparable to the number 
of stars in that galaxy. Close white dwarf-red dwarf binaries are only a small fraction of all 
stars in the Galaxy. If this is also true in M31, then each nova must erupt ~ 4 times, where 
f is the ratio (number of nova binaries/total number of stars). A precise value of f is not 
known, but it is probably in the range 10~ 3 < / < 10 - ®. 

The finding of all the numerical studies cited above, that M $ lO _ ®-lO" lo Jlf0yr~ 1 
is crucial to permit nova recurrence, is central to the hibernation scenario of cataclysmic 
binaries described in §V. 

II. t DIFFUSION DURING PRE-NOVA ENVELOPE ACCRETION 

Many theoretical studies have shown that CNO enrichment of the accreted hydrogen-rich 
envelope is required to produce a fast nova (Starrfield, Truran and Sparks 1978; Fujimoto 
1982). There is, however, no observed correlation between CNO abundances in the ejecta and 
the nova speed class. Shara, Prialnik and Shaviv (1980) found that at least five important 
physical parameters affect the speed class of a nova. These are the white dwarf mass M W( {, 
the hydrogen-rich envelope mass M ent , the envelope metallicity Z en v, the mass accretion rate 
onto the white dwarf M, and the underlying white dwarf luminosity L W( |. These workers and 
MacDonald (1983) suggested that novae of different speed classes lie at different points in 
the (Z env , M env )-plane. More recent work has shown that neither Z C n* nor M env is a free 
parameter that may be chosen arbitrarily, or even independently. Both are determined by 
the accretion phase preceeding the outburst. In current thinking, the fundamental 
parameters which determine the speed class of a nova are M w i,L w fa and M. 
Novae begin to fade (in the optical part of the spectrum) only after they’ve ejected most 
of their envelopes (Bath and Shaviv 1976). M w d, L w j and M determine the white dwarf 
envelope mass at the time of eruption, and the strength of the eruption. This determines 
how quickly the envelope is ejected, and hence how quickly the optical nova fades (i.e., the 
speed class). 

M w 4 , M and L w d change as the cataclysmic binary evolves. Thus the character and 
frequency of nova outbursts should change in a given system over thousands of outbursts 
spanning many millenia. A self-consistent, nova binary evolution simulation over many 
outbursts has not yet been attempted (but see §111.2 for the evolution through two outbursts). 

Prialnik and Kovetz (1984) and Kovetz and Prialnik (1985) have proposed a simple 
mechanism by which CNO enhancement of the accreting hydrogen-rich envelope is achieved. 
This mechanism is diffusion of hydrogen into the outer core of the C-0 white dwarf and 
C-0 into the hydrogen-rich envelope from the white dwarf core. Ignition first occurs in a 
very hydrogen-poor, but not hydrogen-free region in the outer white dwarf core. The large 
amount of energy liberated as soon as ignition occurs causes convection to mix all the matter 
from the ignition zone (Z £ 1) up to the surface (Z ~ Z re i dwarf ~ Z©). 
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The profiles of heavy elements at the onset of TNRs for various models are extremely in- 
structive. Figure 1 (from Kovetz and Prialnik 1985) shows that the sharp boundary between 
hydrogen-rich envelope and hydrogen-absent core is smeared out by the diffusion process. In 
all diffusion models, traces of hydrogen diffuse into en outer core layer comparable in mass 
to the accreted hydrogen-rich envelope. This results in CNO abundances in the convective 
envelope at the onset of the simulated nova outburst from Z = 0.08 to 0.40. 

The diffusion calculations carried out by Kovetz and Prialnik (1985) coupled the diffu- 
sion equations, which include concentration, pressure, and thermal terms to the evolution 
equations for each timestep and mass shell. By using a coarse model grid in M ,M wc i, and L 
Kovetz and Prialnik (1985) were able to determine roughly the simultaneous dependence of 
M env and Z en « on these parameters, as shown in Figure 2. 

The white dwarf mass M W( { was found to be the single most important parameter in 
determining the accreted envelope mass A mace- The depth reached by hydrogen diffusing 
into the core is proportional to the accretion interval Atacc = Arrive M . Thus the total 
amount of mass A mg for which X > 0.005 is also strongly correlated with The ratio 

Amg/ Amacc, which directly determines Z env is almost independent of M W( j, but is strongly 
dependent on M. 

Z env is mostly determined by the accretion rate M in the sense that higher Z env values 
result from lower accretion rates. Low M permits hydrogen to diffuse for a longer time, and 
hence to a greater depth into the white dwarf core. 

The accreted envelope mass A m^c increases with decreasing L tt ,j, because low luminos- 
ity corresponds to low envelope temperature and correspondingly longer time to reach TNR. 
Contrary to previous simulations, Kovetz and Prialnik (1985) found that Arriacc decreases 
with decreasing M. This >s because diffusion provides enough enrichment for a strong TNR 
in hotter regions than were possible without diffusion. Thus very fast novae should occur for 
a combination of high Z env (low M) and low M env . 

The diffusion scenario predicts that all novae show some degree of CNO enhancement as 
a result of diffusion. It also explains, in a natural way, how CNO nuclei can be enhanced by 
very different factors in novae of different speed classes, as shown in Figure 2. 

In summary, white dwarfs of different masses and luminosities, undergoing accretion at 
different rates require very different time intervals to reach TNR. Diffusion can produce the 
wide variety of observed Z env in a wide range of M en t>. It is M w ,^, and M that determine 
Z env and M env at the beginning of a TNR, and hence the speed class of an erupting nova. 

II. S SHEAR INSTABILITY MIXING OF ACCRETING MATERIAL 
POSSESSING ANGULAR MOMENTUM 

The hydrogen-rich material accreted by a white dwarf in a cataclysmic binary possesses 
very high specific angular momentum. There must exist an interface region between the 
outer part of the white dwarf core and the inner part of the accretion disk where material 
is transferred from the disk to the core. Important and difficult questions concerning this 
interface are: 

1) How effective is shear mixing at transporting accreting material into the white dwarf 
core? 

2) Are there processes that can transport matter away from the equatorial region of the 
white dwarf on a timescale short in comparison with that required to build up an envelope 
sufficiently massive to generate a TNR? 
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3) What, if any, are the effects of rotational energy dissipation on the thermal equilibrium 
of the white dwarf envelope, and of the resulting TNR? 

4) Are there any observed or predicted consequences of shear mixing of accreted and 
core material? 

Durisen’s (1977) pioneering study of the accretion of a rotating fluid by a be ry tropic 
model of a white dwarf outlined many of the important physical principles and effects that 
must be considered. In particular, the effects of degenerate electron viscosity and various 
fluid instabilities that could lead to turbulence and enhanced transport of angular momentum 
were considered. Two particularly promising mechanisms outlined were the formation of an 
Ekman layer, and oscillating, non-equilibrium meridional circulations. The Ekman spin-up 
theory is particularly attractive and has been treated in some detail by Livio and Truran 
(1987) and is discussed at the end of this section. 

Durisen (1977) concluded that the overall structure of an accreted envelope depends 
critically on the turbulent redistribution of angular momentum. A schematic mo^jl for the 
white dwarf in a cataclysmic binary is adopted from Durisen (1977) and shown in Figure 3. 
The star consists of a degenerate core in slow, nearly rigid rotation; a mostly degenerate, 
relatively narrow, turbulent, high-shear region; and a nearly Keplerian envelope which may 
or may not have a substantial degenerate subregion. 

Dunsen (1977) suggested that the nonsphericity of old nova shells (McLaughlin 1960; 
Warner 1972; Hutchings 1972) could be due to TNRs in nonspherical envelopes of white 
dwarfs. This remains a distinct possibility, though interaction of the expanding shell with 
the accretion disk and/or the binary companion, magnetic fields, or localized TNRs are also 
possible. Durisen (1977) also pointed out that the high-shear region could provide core- 
envelope mixing and the CNO enhancements required by the Starrfield, Sparks, and Truran 
(1974) nova models. 

Kippenhahn and Thomas (1978) (KT) claimed that accretion belts do form on accreting 
white dwarfs. They stated that the Richardson criterion, in which the effects of temperature, 
angular velocity and molecular-weight gradients are taken into account, is the condition that 
determines the onset of shear-mixing. 

In the presence of a gravitational field, turbulence is hindered by a density distribution 
increasing towards the center of gravity. To maintain turbulence caused by shear motion, the 
kinetic energy of turbulent elements must be larger than the work done by those elements 
against gravity and centrifugal force. This is quantified by the Richardson number R t - (for 
details see Prandtl 1942): R» = work done against buoyant and centrifugal forces/kinetic 
energy of turbulence. Marginal stability occurs in the case Rj = 0.25. KT noted that if 
the timescale of accretion is long compared to the timescale during which the instability can 
mix material, then the involved layer is nearly marginal. Freshly accreted material causes 
further exchange of angular me ntum and of hydrogen between the disk and the white 
dwarf. This makes the last accreted layer only slightly unstable. The bottom of the accreted 
layer moves deeper into the white dwarf and nuclear burning eventually starts. 

An important simplification used by KT was that the condition for marginal stability 
was satisfied separately in the radial and horizontal directions. Furthermore, the existence of 
the belt was assumed not to perturb the structure of the white dwarf. The main numerical 
result of Kippenhahn and Thomas (1978) was that the accretion belt was found to extend 
to latitudes of ±21°, roughly independent of depth. 

MacDonald (19836) studied the effects of instabilities due to differential rotation and 
chemical inhomogeneity in the envelopes of accreting stars. By considering growth rates of 
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nonaxisymmetric perturbations MacDonald (19836) claimed that freshly accreted material 
rapidly suffers dynamical instabilities and covers the entire white dwarf. Mixing inwards, 
however, towards the white dwarf core occurs on a much longer (thermal) timescale. 

The ratio of soft X-ray flux to visual flux for many old novae is several orders of mag- 
nitude less than predicted by Pringle’s (1977) study of the boundary layer between the 
accretion disk and white dwarf (Ferland et al. 1982). Perhaps the simplest solution to this 
problem (suggested by Ferland et al. 1982) is that of a boundary layer comparable to the 
size of the white dwarf rather than a thin belt at the white dwarf’s equator. A larger sur- 
face area permits a lower effective temperature in the boundary layer for the dissipation 
of rotational energy. Unobservable EUV radiation, rather than soft X-rays, would then be 
emitted. MacDonald (19836) suggested that shear-induced turbulence spreads the thermal 
energy (produced by the dissipation of rotational energy) evenly over surfaces of constant 
effective gravity on a dynamical timescale. “The missing boundary layer” becomes, in this 
interpretation, evidence that shear-induced turbulence mixes accreted matter on a dynamical 
timescale over the entire white dwarf. 

Sparks and Kutter (1987) and Kutter and Sparks (1987) carried out extremely detailed 
numerical simulations of accretion and thermonuclear runaways on white dwarfs. They too 
assumed that the angular momentum distribution of accreted material on a white dwarf was 
in marginal stability and was determined by the Richardson criterion. The distribution of 
mass they considered had a half width of 27.5° at half intensity, covering 42.3% of the white 
dwarf’s surface. This is more of a sash than a belt. In all cases studied, Sparks and Kutter 
(1987) failed to obtain a sufficiently strong nuclear runaway for nova-like mass ejection. In 
their radial accretion models gravitational compression heating triggered the TNR too soon. 
When they included the effects of angular momentum on the accreted matter, the support 
due to the centrifigal force lowered the mechanical confinement to the burning regions too 
much to permit a strong TNR and mass ejection. 

Sparks and Kutter (1987) considered several possibilities to increase the pressure at the 
base of the accreted material during TNRs. One suggestion was that only white dwarfs 
significantly more massive than 1 Mq can generate classical novae. However, the estimated 
masses of white dwarfs in nova binaries is generally 1 Mq or less (Warner 1976, Ritter 
and Burkert 1986; but see Livio and Soker 1984). Another possibility is that the angular 
momentum of the accreted matter might be removed by tidal torques’ braking action exerted 
by the disk or by the secondary star, or possibly by mechanical coupling between the accreted 
matter and the disk. This could permit mass to penetrate the core without carrying in 
angular momentum. The third possibility is the interaction (during a TNR) of convection 
with layers of material whose horizontal velocities have a gradient. 

Livio and Truran (1987) considered the generation of a very thin turbulent layer on a 
white dwarf surface dynamical time scale (seconds) as soon as the accretion disk first comes 
into contact with the white dwarf. The concluded that on a turbulent, viscous tim< scale an 
almost rigidly rotating thin layer spreads over the entire white dwarf surface. As the thin 
rotational shear layer (the “Ekman layer”) becomes turbulent it produces another shear layer, 
with the process repeating itself. This mechanism transfers angular momentum towards the 
poles and to the body of the white dwarf much more efficiently than the Kippenhahn-Thomas 
model. Spreading the angular momentum over the white dwarf’s surface and into the core 
reduces the centrifigal force in the outer layers and premits nova-like TNRs in the Sparks 
and Kutter scenario. 
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Researchers in the field are still debating whether diffusion or shear mixing is the domi- 
nant mechanism in producing CNO and heavy element enrichments in nova ejecta. A simple 
and elegant test has been proposed by Livio and Truran (1987). The test is based on the 
conclusion of Webbink et al. (1987) that the outbursts of the recurrent novae T Pyx and U 
Sco (and by implication V394 CrA) result from TNRs on the surfaces of very massive white 
dwarfs. Because of the short recurrence times of these objects (1-2 decades) diffusion (which 
requires centuries for significant mixing) cannot be operative in them. Any-heavy element 
abundances seen would have to be due to shear mixing. Conversely, a complete abs'*' 
of heavy elements in the ejecta spectra of all three recurrent novae would imply tha* . " 
mixing is rather unimportant on all nova white dwarfs. Diffusion would remain th 
contender as the enrichment mechanism in non-recurrent novae. 

T Pyx is due for an eruption in the next few years and should be carefully momt.-; 
spectroscopically for heavy element enrichments. U Sco is not enriched in CNO abundances 
relative to the Sun (Williams et al. 1981). However, it may be heavily enriched in helium. 
Shear mixing might have occurred with a helium-rich layer left behind from the previous 
outburst. Th* observed high ejection velocities remain unexplained. 

A photometric and spectrographic study of the second known eruption of V394 CrA has 
been reported by Sekiguchi (1988). The amplitude of the outburst (AV « 11 magnitudes) 
is comparable to that of classical nova. The decline rate (3 magnitudes in 5 days), high 
ejection velocities and spectra resemble those of U Sco. 

II.4 THE EFFECTS OF HEATING FROM THE BOUNDARY LA YER 
ON ACCRETION MODELS FOR NOVAE 

Once matter in the accretion disk of a cataclysmic binary reaches the environs of the 
white dwarf its dynamical state must change quickly. While in the accretion disk the material 
orbits the white dwarf with Keplerian velocity. There must be strong thermal and dynamical 
interaction between the outer part of the white dwarf and inner part of the accretion disk in 
order for accreting material to transfer from the latter to the former. The structure of this 
boundary layer has not yet been solved. However, Shaviv and Starrfield (1987) have made 
the point that boundary layer heating of the outer white dwarf layers must be important, 
and must be included in accretion models of cataclysmics. 

The virial theorem implies that at least half the boundary layer luminosity of material 
accreting onto a white dwarf is radiated, while up to half is transported into the surface layers 
of the white dw'arf by the accreted material. The effects of this accreted energy significantly 
affect the evolution to a TNR and nova outburst. In particular, Shaviv and Starrfield (1987) 
found that inclusion of boundary layer energy makes it more difficult to produce a nova 
outburst for a given set of parameters (M, M w ^i e dwarf % ^white draw f) usually chosen in 
studying novae. This is because boundary layer heating was found to produce a nuclear 
burning region extending through all or much of the accreted envelope. This burning forced 
the entire accreted envelope to become convective. Because the degeneracy of material in 
the envelope was significantly lowered by the heating, the strengths of TNRs were also 
dramatically reduced. 

Shaviv and Starrfield (1987) found their results to be strongly dependent on the degree 
of boundary layer heating and on the ratio of mixing length to scale height. In the words of 
the authors... “we have complicated an already cloudy situation”. Nevertheless, it is clear 
that heating from the boundary layer must be included in realistic accretion models. 
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The reader should not be mislead into thinking that TNRs are impossible to produce 
when the effects of boundary layer heating are included. We note that both the 1.35 A/© 
accreting models of Starrfield, Sparks and Shaviv (1988) and the l.lAf© iron white dwa'-f 
models of Prialnik, Kovetz, and Shara (1988) produce strong TNRs even when boundary 
layer heating was included. The latter authors concluded that the maximum M that could 
be tolerated, and still produce a strong TNR on a white dwarf of given mass, had to be 
decreasedby about an order of magnitude because of boundary layer heating. Thus accretion 
rates £ 10 ' 10 Af©yr _1 are theoretically necessary for TNRs when boundary layer heating 
is considered. 


III. THE OUTBURST STAGE 

III. I LOCALIZED THERMONUCLEAR RUNAWAYS (LTNR) 

Most studies of thermonuclear runaways on the surfaces f white dwarf stars have made 
the simplifying assumption that the runaways are completely spherically symmetric. This 
is a reasonable first approximation in trying to understand what the dominant physical 
processes are in the runaway. However, even simple estimates (described below) show that 
TNRs usually start at a point somewhere in the surface layers of an accreting white dwarf 
star. How likely it is for the thermonuclear runaway to spread from the point, and how 
quickly spreading can occur are much more involved problems. 

The first detailed study of localized thermonuclear runaways was conducted by Mitro- 
fanov (1980), who was concerned with TNRs on magnetic white dwarfs. Mitrofanov sug- 
gested that the outbursts of dwarf novae result from the unstable thermonuclear burning of 
hydrogen on the surfaces of magnetic degenerate dwarfs. He noted that with a magnetic field 
of 10«- 8 gauss a white dwarf would be able to confine hydrogen-’‘ich plasma to accretion 
columns. A temperature inhomogeneity on an accreting white dwarf’s surface is then gen- 
erated as accretion heats only a small part, of the surface. Furthermore, a strong magnetic 
field leads to heterogeneities in the heat conductivity. A typical area exposed to accretion 
was crudely estimated to be about 10 14 " 15 cm 2 , or a few percent of the surface area of the 
white dwarf. 

Mitrofanov pointed out that a competition takes place between the rate of local accumu- 
lation of material at a magnetic pole and the rate at which material spreads away from the 
pole over the surface of the white dwarf. If local accumulation is rapid then ignition might 
occur in a small area on the surface of the white dwarf. He estimated the miss of burning 
hydrogen to be approximately 10 23 to 10 24 grams. Thus if 10 18 ergs gm _1 are released dur- 
ing the TNR (an amount of energy per gram similar to that released in a nova eruption) 
then the total dwarf nova outburst energy is about 10 39 to 10 49 ergs. This agrees with the 
energetics of dwarf nova outbursts. 

The recent demonstration that Nova Cygni 1975 is a strongly magnetic white dwarf (see 
§VI.2) has cast serious doubt on the picture of Mitrofanov. The magnetic field of Nova Cygni 
1975 did not confine most of the hydrogen to small areas at the polar cap of the star. The 
total energetics of Nova Cygni demand that about 10 -S M© of hydrogen was accumulated 
on the white-dwarf surface without being burnt to helium (10~ 8 M© at a time) in successive 
dwarf nova eruptions. While Mitrofanov’s hypothesis of dwarf nova eruptions being caused 
by thermonuclear runaways is probably incorrect, his paper was the first to point out that 
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localised thermonuclear runaways should occur on white dwarf stars under certain special 
circumstances. 

The next studies of localized runaways were carried out by Shara (1982) and by Fryx- 
ell and Woosley (1982). Shara (1982) considered the energy transport inside a degenerate 
hydrogen rich envelope on a white dwarf. lie demonstrated that only when thermal conduc- 
tivity is the dominant mode of energy transport inrde a hydrogen rich envelope can a local 
TNR occur. The thermalization time scale between adjacent mass elements AR in size is 


tj<H ~ 10 4 


[ATT 

2 

L10 7 . 



1 -(- 


0.45 (T/10 7 ) 2 (l0 4 /p) 2 ] 1 ) years, 


where T is the temperature, p is the density and If is the Fermi temperature. For typical 
white dwarf envelope values (p ~ 10 4 gm cm -1 and T ~ 15-20 x 10® °K) a localized TNR 
is able to occur long before the envelope can reach thermal equilibrium. 

A VNR should start in the initially hottest part of the envelope covering only a small 
part of the white dwarf surface area. In some cases neither conductive, convective, radiative, 
advective, nor acoustic transport of energy from the site of the LTNR seem capable of heating 
all or much of the surrounding hydrogen-rich envelope on a time scale comparable with the 
eruptive phase of the TNR. 

Temperature inhomogeneities are crucial in making a localized rather than a global TNR. 
Four mechanisms for setting up temperature gradients were suggested by Shara (1982). 
These are the degenerate dwarf’s magnetic field, rotation, non-spherically symmetric accre- 
tion of matter, and the red-dwarf tidal force. 

As discussed in §11.3 the accretion of matter by the white dwarf is certainly rot spherically 
symmetric. Matter is accreted via an accretion disk, nence the deposition of matter and 
turbulent energy on the white dwarf outer layers occurs in a toroidal region around the 
equator of the degenerate dwarf. More compressional (Pd v ) work is done on material in this 
toroidal region than on material at higher latitudes on the white dwarf surface. Temperature 
differences of a few percent or even larger over regions 10 7 8 cm apart are straightforward 
to produce. Hydrogen-rich envelopes must be massive for local TNRs to occur because only 
then does thermal conductivity dominate radiative energy transport. 

If a localized runaway does begin in the envelope of a degenerate dwarf, what are the 
logical consequences? In particular, will a localized runaway remain localized, or will it 
spread in a short time to engulf the entire dwarf envelope? 

Once a TNR (local or global) starts on a white dwarf, radial transport of energy is 
convective. Determining the method of meridional transport of energy is much more difficult. 
Becaus' the dominant energy transport mechanism is not certain, neither is the timescale of 
energy transport. As a result, it’s not certain how long localized runaways remain localized. 
In the absence of turbulent mixing or horizontal convection, Shara (1982) claimed that it 
would take years for a diffusively propagated burning wave to cover the surface of a white 
dwarf. 

A localized TNR has a built in “safety valve” not available to a spherically symmetric 
TNR. In the latter case, the only direction in which matter can expand is radially. In the 
localized case, horizontal expansion is also possible. In particular, as material expands quasi- 
hydrostatically out of a cylindrically shaped LTNR, that material must spread horizontally, 
away from the site of the localized runaway. This is sketched in Figure 4. A volcanic out- 
pouring of hot, hydrogen-helium “lava” onto the white dwarf surface h the obvious analogy. 
Convective entrainment and turb’ilence were not included in these qualitative scenarios, and 
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any estimate of the spreading time scale would depend on these effects. Fryxell and Woosley 
(1982) suggested that these processes spread the runaway over the entire dwarf surface in 
days or weeks. 

Orio and Shaviv (1988) have recently extended the analytic estimates made by Shara 
(1982) and by Fryxell and Woosley (1982). They included the accretion process of the 
hydrogen-rich material onto the white dwarf and they extended the calculations to areas 
comparable to the white dwarf surface. The most importantd physics introduced by Orio 
and Shaviv (1988) is consideration of the fraction of the radiative flux tj distributed in the 
meridional direction. By taking into account the curvature of the white dwarf’s surface, Orio 
and Shaviv (1988) also allowed for temperature inhomogeneities on scales of the order of the 
ra ’lus of the whit~ dwarf. 

The general trends found confirm many of Shara’s (1982) results. In particular, the 
result that TNRs which start at a point can remain local is confirmed and extended by Orio 
and Shaviv’s (1988) work. They find that for a minimum temperature differential AT a very 
small TNR will spread into a somewhat larger one, but not quickly cover the entire white 
dwarf surface. 

Numerical simulations indicate that LfNRs are most likely to occur on white dwarfs 
accreting at high rates, on hot massive white dwarfs, or on cool intermediate mass white 
dwarfs. Meridional transport of energy lengthens the propagation time of a local TNR. 
A months or years-long brightening before the rapid, dramatic rise to maximum light is 
frequently observed in classical novae. Orio and Shaviv (1988) suggest thai this rise to 
maximum light might be due to meridional transport of energy slowing down the TNR. The 
pre-nova EUV radiation from the white dwarf could cause weeks or months of heating and 
brightening of the red dwarf (seen in visible light). 

Orio and Shaviv (1988) have demonstrated how complex LTNRs can be; until convection 
and turbulence are included in realistic, three dimensional simulations of local TNRs, the 
effects and behavior of LTNRs on white dwarfs will remain rather poorly understood. 

III. 2 THE EVOLUTION OF A CLASSICAL NOVA MODEL THROUGH 

A COMPLETE CYCLE 

The simulated evolution of even a single thermonuclear runaway on a white dwarf surface 
can take hours of time on a large computer. The coupled equations of hydrodynamics and 
nuclear energy generation and transport, diffusion, and mass accretion make any attempts 
at accurate analytic solutions hopeless. 

Until very recently, the complexities of the numerical simulations have prevented even 
the simulation of two successive runaways on the same white dwarf. Prialnik (1986) has 
now followed the evolution of a classical nova through a complete cycle. The accretion, 
outburst, mass loss, decline, and resumed accretion, leading to a second outburst, were all 
calculated with a hydrodynamic code, which included diffusion and an extensive nuclear 
reaction network. The model’s luminosity history is shown in Figure 5, and that of the 
envelope is shown in Figure 6. The white dwarf chosen was a 1.25 M© C-0 star, while the 
accretion rate was taken to be M =■ lO' ^Afoyr -1 . Matter with close to solar composition 
was accreted onto the white dwarf. The exact parameters of the model are less important 
than the result of the simulation which is that the second accretion phase resembled the 
first one very closely, and the second outburst was almost identical to the first. While this 
has been suspected by many workers in the field for many years, this work was the first 
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clear demonstration of the correctness of this statement. The implications for the long term 
evolution of classical novae, and of successive nova eruptions, are extremely important. 

Many of the results of Prialnik (1986) of the first stage of accretion through outburst 
and decline have already been seen in previous TNR simulations. The low accretion rate 
chosen permitted significant amounts of diffusion of hydrogen-rich material into the C-0 
white dwarf and the resulting enrichment over solar (by an order of magnitude) in Z lead to 
a fast nova. The composition of the ejecta and the light curve and ejection velocities were 
in good agreement with e.g. Nova Cyg 1975 and Nova Cyg 1978. Within a few centuries, 
the white dwarf cooled to its preoutburst structure, with one very significant difference. A 
low-mass but helium-rich envelope was left on the white dwarf. It was this layer which raised 
the helium abundance of the ejected matter of the next nova outburst (see the next section). 

Because of diffusion and the resulting enrichment in CNO nuclei of the accreted hydrogen- 
rich envelope, the white dwarf mass was found to decrease slightly. In addition, the red dwarf 
mass decreased secularly because of mass transfer to the white dwarf. Thus the orbital 
period, which depends only on the secondary mass (Patterson 1984) should decrease. In 
the event of very low mass-rate accretion, like that considered in Prialnik’s work, the white 
dwarf ’8 intrinsic luminosity decreases steadily and the occurrence of flashes does not prevent 
the normal cooling of the white dwarf. However, for higher mass accretion rates, where 
heating from the boundary layer must be included, the opposite may be true. Whether 
the white dwarf luminosity increases or decreases is critical to understanding the long term 
evolution of a nova binary. As Prialnik (1986) notes “a significant decrease in the white 
dwarf temperature would result in an appreciable increase of the accretion time prior to 
outburst and, consequently, larger envelopes and ejected masses.” The opposite statement 
is also true. Thus, despite this important first step in understanding the long term evolution 
of classical novae, a simulation of several successive outbursts which include a self-consistent 
determination of the mass transfer rate throughout the outburst cycles is still needed to 
determine the long-term trends in nova eruption behaviour. 

III.S THE HELIUM/HYDROGEN RATIO IN THE EJECTA OF 

CLASSICAL NOVAE 

Highly enhanced CNO abundances are not always required to produce a TNR and a 
classical nova (Prialnik, Shara, Shaviv 1978, 1979; Sparks, Starrfleld and Truran 1978; Nariai, 
Nomoto and Sugimoto 1980). Observations of the spectra of novae ejecta, and many TNR 
simulations have shown that a variety of nova speed classes can be produced by a wide variety 
of CNO isotope abundances. We note that various combinations of large overabundances, 
relative to the Sun, of helium, carbon, nitrogen, oxygen, neon, magnesium and perhaps iron 
have all been seen in various novae (see Truran and Livio 1986 for a recent summary). 

There seem to be only three reasonable ways to create this overabundance (Prialnik 
1988). These are: 

1) the accreted matter must itself be extremely rich in helium, CNO, and perhaps heavier 
isotopes; or 

2) during the outburst, the heavy elements are created by nucleosynthesis; or 

3) the accreted and core materials are mixed by diffusion and/or shear mixing during the 
accretion process. 

The mass transferring red dwarfs in cataclysmic binaries in the solar neighborhood should 
be composed of solar-like material. Furthermore, no indication of abundance pecularities in 
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the accretion disk spectra of nova-like variables (widely believed to be nova progenitors) has 
ever been reported. Thus the first possibility can be discounted. Nucleosynthesis of large 
amounts of heavy isotopes can also be discounted. Less than 1% of the hydrogen envelope 
need be burnt in order to provide the energy of a nova outburst. 

If mixing is the only viable explanation for high metallicity and the violent mass ejection 
seen in novae, and heavy isotope production during thermonuclear burning is not important, 
then the ratio of He/H should be that of the accreted material, t.e., solar. The observations, 
recently summarized by Prialnik (1988) show that enhanced He/H ratios are very common. 
Three outstanding examples are RR Pic (nova 1925), DK Lac (nova 1950), and HR Del 
(nova 1967). In each of these three objects the ratio Y/X is approximately 1.0, while the 
. ~>lar value is 0.34. An inescapable conclusion seems to be that the outer layers of the white 
dwarf are very rich in helium. 

A significant complication to this simple picture, however, is that large CNO excesses 
usually accompany the helium excesses. Thus mixing must be able to bring hydrogen- rich 
material down through the top helium layer into the white dwarf C-0 core. The helium 
layer thickness must, therefore, be comparable to the hydrogen-rich envelope mass, i.e. of 
order 10 -5 -10“ 4 M©. A thin primordial shell of helium would be quickly eroded by just a 
few nova outbursts, and would almost never be seen as enriched helium in the ejecta of most 
novae. 

An elegant alternative mechanism, proposed by Prialnik (1988) is that a thin helium-rich 
mantle forms just after each nova outburst phase. This occurs because only a fraction of 
the hydrogen-burning envelope is ejected during each nova eruption. Some fraction of the 
mass, 30%, is left behind and much of this material burns to helium. Such behaviour 
has been demonstrated in TNR simulations which eject mass. The mass accreted before 
the next eruption constitutes between 40-80% of the total envelope mass. Diffusion, and/or 
shear mixing enrich the envelope both in helium (left over by the last eruption) and in C-0 
or O-Ne-Mg from the white dwarf. 

In the case of a C-0 white dwarf, most of the carbon dredged up from the core is burned 
into nitrogen by the eruption. The TNR temperatures are generally too low for oxygen to 
be affected (Starrfield et al. 1978). The C/O and high He/H ratios observed in novae ejecta 
provide strong support of tin scenario of the depletion and regeneration of a thin helium 
shell during each nova cycle. 

III. 4 THERMONUCLEAR R UNA WA YS ON THE S URFA CES OF 

OXYGEN-NEON-MAGNESIUM WHITE DWARF 

As long as the temperature in the L/drogen-rich envelope of a white dwarf, which is 
undergoing a thermonuclear runaway, does not exceed 4 x 10® °K “breakout” from the CNO 
cycle does not occur (Truran 1982; Wiescher et al. 1986). The C, N, and O isotopes present 
act only as catalysts; they are neither created nor destroyed in the CNO cycle burning. Only 
when T > 4 x 10® °K can the Ne-Mg-Al-Si group of nuclei be created (i.e., breakout). No 
single published TNR model to date, including those dealing with white dwarfs approaching 
the Chandrasekhar mass, reaches or exceeds 4 x 10® °K. Thus, any element heavier than 
oxygen, seen in nova ejecta, must either have been transferred into the white dwarf envelope 
from the red dwarf companion, or somehow have been dredged up from the underlying white 
dwarf. 
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One of the first quantitative reports of high element abundance was that by Ferland and 
Shields (1978) of a twenty-fold enhancement of the abundance of neon relative to solar in 
the ejecta of nova Cygni 1975. Since then, even more dramatic enhancements have been 
reported by Snijders et al. (1984) for nova V1370 AQL 1982, and by Williams et al. (1985) 
for nova V693 CrA 1981 (see Figure 7). Andrillat and Houziaux (1985) observed very strong 
Mg II emission both in IUE and optical spectra, for the nova Vul 2 1984. 

Ferland and Shields (1978), Law and Ritter (1983), and Williams et al. (1985) have 
all suggested that the high neon abundances seen in these novae imply the existance of 
underlying O-Ne-Mg white dwarfs. Just as the high abundances of C, N, and 0 in many 
nova ejecta imply the ablation of the underlying white dwarf, so does the presence of oxygen- 
neon-magnesium enriched material imply the erosion by the nova process of the underlying 
white dwarf. 

Stellar evolution calculations predict that stars with initial masses in the range 1-8 Af© 
will evolve to yield C-0 white dwarfs. Those of initial mass 8-12 Mq will evolve further to 
yield ONeMg white dwarfs (Barkat, Reiss and Rakavy 1974). A simple Salpeter (1955) initial 
mass spectrum predicts the eventual formation of 35 C-0 white dwarfs for every ONeMg 
white dwarf. Iben and Tutukov (1984, 1985) arrived at a similar ratio, after including possible 
effects of binary evolution and mass transfer onto the white dwarf. 

Because at least four of the dozen recent brightest novae have exhibited strong neon 
and/or magnesium enhancements, the frequency of oxygen-neon-magnesium white dwarfs 
appears significantly larger than that predicted by the above theoretical considerations. 
However, Truran and Livio (1936) have shown that selection effects are very important in 
determining which kinds of white dwarfs are seen in nova eruptions. ONeMg white dwarfs 
may be more massive on average than C-O white dwarfs, and, therefore, are predicted to 
erupt with lower mass envelopes. If this assumption is correct then ONeMg white dwarfs 
accumulate critical mass envelopes and suffer nova eruptions at shorter recurrence intervals. 
Thus, despite the fact that ONeMg white dwarfs occur in only ~ 3% of all cataclysmic 
binaries, their large masses ensure that they account for ~ 1/3 of all observed classical 
novae. 

The first hydrodynamic simulations of oxygen-neon-magnesium novae have been carried 
out by Starrfield, Sparks and Truran (1986). These authors studied thermonuclear runaways 
on the surface of a 1.25 M© white dwarf. Accretion rates of 10 _8 -10~ 10 M©2/r -1 , as well 
as different chemical compositions in the accreted matter were tested. Mixtures with O 16 
abundances of 0.25 and 0.50 were used to simulate accretion onto oxygen-neon-magnesium 
white dwarfs. The results of these simulations were extremely violent TNRs and nova models 
which exceeded the Eddington luminosity L ^dd maximum. Ejection velocities of over 
4 x 10 8 km s -1 were achieved. While the ejection velocities are higher than those in C-0 
white dwarf models, the ejected masses are significantly lower. 

More realistic accretion scenarios, including diffusion and shear mixing, should be at- 
tempted to try to answer some of the remaining difficulties with the oxygen-neon-magnesium 
white dwarf scenario. In particular, is it possible to make a nova eject roughly comparable 
fractions of hydrogen, helium, CNO and neon as seen in V693 CrA (nova 1981)? At least 
as puzzling is the appearance of significant enhancements of iron relative to solar in several 
recent novae. As t: Larrfield (1987) has emphasized, high iron abundances are not predicted 
by TNR scenarios on either C-0 or ONeMg white dwarfs. A first attempt at explaining these 
high-iron observations has been made by Prialnik, Kovetz and Shara (1988). The results are 
described in the next section. 
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III. 5 THERMONUCLEAR RUNAWAYS ON THE SURFACE <? OF 
IRON WHITE DWARFS 

As described above, the large enrichments seen in several novae of oxygen, neon and 
magnesium, relative to solar, have suggested a population of ONeMg white dwarfs and 
cataclysmic binaries. Significant overabundances of iron have also been seen (in N Muscae 
1983 and N Aquilae 1982). Conventional thermonuclear runaway theory does not predict 
the formation of la: 0 e amounts of iron during nova eruptions. This has prompted Prialnik, 
Kovetz, and Shara (1988) to study accretion onto, and nova events on the surfaces of Fe 
white dwarfs. Half the free-fall energy of the accreted matter was assumed to be transmitted 
to the accreting star via the accretion disk boundary layer (a heating process ignored in most 
previous simulations). Diffusion between H, He, CNO, and Fe was taken into account. 

The material which was accreted onto the iron white dwarf was assumed to be of solar 
composition. During the accretion phase hydrogen was found to diffuse deeper into the Fe 
core than its CNO catalysts. Thus a TNR did not begin at the deepest levels penetrated 
by the hydrogen, but further out where the CNO abundance was highest. This natural 
enhancement of the CNO abundances does serve to trigger TNRs in every case studied, but 
these overabundances were smeared out by convection as soon as the TNR began. Further- 
more, accretion rates M of 10 — 10 — 1 ®, and 10 -11 Mo/yr produced qualitatively different 
runaways because of different envelope masses and thermal structures at the time of TNR. 

Large iron enhancements (over solar) were obtained, particularly at M = 10 _1 ® M©yr -1 . 
Generally, a high degree of contamination of nova ejecta by core material is obtained either 
from rapid accretion onto a non C-0 white dwarf or from slow accretion onto a C-0 white 
dwarf (if diffusion followed by convection is the mixing mechanism). In the former case, this 
is because the admixed core material does not participate in the thermonuclear burning, 
and because hydrogen will only diffuse to a Lagrangean mass depth in the white dwarf core 
of a few times the mass of the hydrogen-rich accreted shell. A long period of slow mass 
accretion yields a more massive envelope, but not one more enriched in iron. In the case of 
a C-0 white dwarf, the admixed carbon and oxygen enhance the TNR so slow accretion is 
preferable to obtain a high degree of C-0 enhancement. 

Mass loss from the iron white dwarf was found to occur in two distinct stages. The 
episodes of mass loss occurred several days or weeks apart as the erupting star expanded, 
ejected matter, contracted and then expanded again, shedding more mass. RR Pic (nova 
1925) showed multiple mass loss episodes, and probably was iron-rich (Payne-Gaposhkin, 
1957). 

The maximum temperature obtained — 2.85 x 10 8 o K — was not high enough for heavy 
elements (CNO, Ne) to be produced. Thus there was no breakout from the CNO cycle. The 
authors concluded that, in order to account for the large iron abundances detected in the 
ejecta of some classical novae, within the framework of conventional TNR theory, the white 
dwarf progenitor should be rich in iron as well as CNO and Ne. 
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///. 6 THERMONUCLEAR R UNA WA YS ON THE SURFA CES OF 

WHITE DWARFS NEAR THE CHANDRASEKHAR LIMIT 

Recurrent novae have been exhaustively reviewed by Webbink, Livio, Truran and Orio 
(1987). These authors concluded that the only two near-certain thermonuclear runaway 
powered recurrent novae were T Pyx and U Sco. The extrordinarily ahort recurrence times 
of these two novae, (eight years in the case of U Sco!) imply that the underlying white dwarfs 
must be more massive than ~ 1.3Af© (Nariai and Nomoto 1979; Fujimoto 1982; MacDonald 
1983) if it really is a TNR that powers these objects. This is because the high surface gravity 
of a massive white dwarf can produce a TNR with a much lower envelope mass (which can 
be accreted quickly) than that on a low mass white dwarf. 

The first simulations of TNRs on the surfaces of white dwarfs near the Chandrasekhar 
limit were carried out by Shara (1981) to see if “breakout” from the CNO cycle could be 
obtained (i.e., if Ne or Mg nuclei could be created during hot hydrogen burning). Tempera- 
tures of 3-3.5 xIO 8 # were the maximum found in these simulations, and “breakout” (which 
requires T ~ 4 x 10 8 °K) did not occur. 

Starrfield, Sparks, and Truran (1985) studied a 1.38Af© white dwarf accreting matter 
with solar composition at 1.7 x 10~ 8 M©j/r -1 and 1.7 x 10 - ®Af©yr -1 . In the former case, 
only ~ 33 years were required to reach the peak of a TNR which ejected 6% of the accreted 
envelope. The remaining material burned to helium in ~ 2 years and settled back onto 
the white dwarf. The light curve, accretion time required to reach outburst, and quantity 
of mass ejected during this simulation agreed well with the observed outburst properties of 
nova U Sco 1979. 

The lower accretion rate simulation took over 10 3 years to reach the burst phase of the 
evolution. Again, only a small fraction (13%) of the accreted envelope was ejected at slow 
velocities. Thus in both cases the mass of the white dwarf was found to grow toward the 
Chandrasekhar limit. 

To try to explain the extraordinary eight year recurrence time of U Sco (eruptions ob- 
served both in 1979 and in 1987), Starrfield, Sparks and Shaviv (1988) simulated the accretion 
of 1.1 x 10 _ ®Af©i/r _1 onto a 1.35M© white dwarf. Incredibly, a TNR was achieved in only 
~ 2.6 years, and it was able to eject 4 x 10 -7 M© moving in excess of 400 km s 1 . The 
same important conclusion as before was reached vis a vis the ultimate fate of the white 
dwarf. Its mass is inexorably increasing, as only a small fraction of the accreted envelope is 
ejected during each eruption. Thus, the authors predict that U Sco will absorb enough mass 
in a rather short length of time (~ 10 5 years) to exceed the Chandrasekhar limit, and will 
erupt as a type I supernova. 



IV. THE EJECTION PHASE 


IV. 1 THE STEADY STATE, CONSTANT LUMINOSITY, CONTINUOUS 
WIND EJECTION MODEL OF NOVAE 

It is a remarkable fact that all erupting novae display the same sequence of spectral 
behavior as the light declines (McLaughlin 1960). While it is certainly true that individ- 
ual novae display unique characteristics in their outburst properties (perhaps due to vary- 
ing binary separation, masses, abundances, etc.) the overall similarities are so strong that 
some universally applicable physical conditions must hold in the ejecta of novae. This basic 
philosophy has been adopted by Bath and Shaviv (1976) in their study of the outflowing 
photosphere of a classical nova. 

All novae near peak brightness display the approximately Eddington luminosity of a ~ 1- 
2 Mq object. This led Bath and Shaviv (1976) to suggest the following set of conditions. 
Radiation pressure in the optically thick continuum of the envelope was assumed to continu- 
ously accelerate matter away from the underlying white dwarf. Electron scattering was taken 
to be the dominant opacity mechanism. The total bolometric luminosities of erupting novae 
were assumed to be constant in time throughout the optical outburst. (UV, infrared and 
X-ray observations have since confirmed this suggestion). Finally, the A-F supergiant-like 
spectra of novae near maximum light (with L~ 10 4 Lq) imply photospheric temperatures 
~ 10, 000 ° K and thus radii of ~ 10* 2 cm. Such radii are ~ ten times larger than the binary 
separation between the white and red dwarfs. Thus frictional dissipation must result from 
binary motion within the common outflowing envelope. The optical decline of a nova was 
found by Bath and Shaviv (1976) to be due to an increasing bolometric correction as the 
wind thinned out, and the photosphere receded to hotter, deeper layers. 

On the basis of their model Bath and Shaviv (1976) were able to determine the mass 
outflow rates from classical novae as M ~ 10 2 *-10 22 <7m a - * at maximum light. Most im- 
portant of all, they showed that the fractional decline from maximum light in the optical 
depends only on the physical conditions in the outflowing photosphere. This accounts for the 
observed strong correlation betwen “decline stage” (t.e., brightness decrease from maximum 
light) and spectral type observed in classical novae. 

Optically thick winds in classical novae were examined in more detail by Bath (1978) and 
by Ruggles and Bath (1979). These studies showed that as the wind thins during the optical 
decline of classical novae, the photospheric location and hence temperatures rises as shown 
in Figure 8. Ions of increasing ionization potential are then photoionized by the increasingly 
hotter emergent spectrum. Bath’s (1978) sequence of models with progressively thinner wind 
and thus higher effective temperatures yielded the decline stage at which oxygen and nitrogen 
could be photoionized by the underlying continuum. The predicted stages at which various 
ionized species are produced by photoionization by the increasing ultraviolet continuum are 
in good agreement with those observed in average novae and are shown in Figure 9. Ruggles 
and Bath (1979) further showed that the outflowing wind must be adiabatic deep within 
the photosphere and in the region near the critical (sonic) point, while small deviations 
from adiabatic flow determine the transsonic structure in the outer envelope. Only flows 
which become supersonic close to the white dwarf surface can produce the ejection velocities 
observed in novae. While radiation pressure is dominant throughout the envelope, most of 
the luminosity of the eruption is at first carried away from the white dwarf by advection. In 
the outer part of the envelope most of the luminosity is transported outwards by radiation. 
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Friedjung (1987) and Taylor et al (1987) have emphasized the importance of the colli- 
sion of different velocity materials in nova ejecta. Friedjung pointed out that the principal 
absorption system, which appears to represent most of the mass ejected by a classical nova, 
could be explained in the following way. Relatively slowly moving material, ejected during 
the premaximum phase of a nova outburst, is overtaken by more rapidly moving ejecta in the 
form of a high velocity wind seem after maximum and which is responsible for the emission 
lines seen in the “diffuse enhanced” and “Orion” spectrum (McLaughlin, 1960). This could 
lead to the production of a relatively dense shell, which Friedjung suggests is the matter 
responsible for the “principal system”. 

Taylor et al (1987) suggest the same model to explain their radio observations of Nova 
Vulpeculae 1984 number 2. The general features of the radio light curve were explained by 
a wind with a mass loss rate M~ 10~ 6 Af©yr -1 , continuing for a year after optical outburst. 
The early radio emission, which has a brightness temperature in excess of 10 5 °K , was 
modeled as a strong shock propagating outward through the principal ejecta of the nova. 
The source of the shock is interaction of the s’ower, earlier ejected material with a later, high 
velocity wind from the white dwarf. 

Detailed models to account for the spectra of nova during outburst have not yet been 
attempted. The wealth of high quality ultraviolet, visible and infrared spectra of novae 
in outburst are just beginning to be qualitatively understood. A detailed quantitative un- 
derstanding and simulation of outburst spectra may require an effort comparable to that 
invested by theorists over the last two decades in simulating thermonuclear-runaways. 
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IV. 2 DUST IN NOVA SYSTEMS 

A comprehensive review of observations and theory of dust in nova systems has been 
given by Bode and Evans (1988). I will discuss only a few of that review’s highlights and 
mention more recent IRAS results. 

The seminal observations were those of FH Serpentis (1970) from 1 to 22 / im by Geisel, 
Kleinmann and Low (1970). These authors observed the visible and infrared fluxes to decline 
together for seven weeks before unexpectedly starting to differ dramatically. The visual light 
curve underwent a steep decline, while the IR flux rose steadily until it dominated the nova 
light output. Similar behavior in Nova Vul 1976, and the subsequent IR decline is shown in 
Figure 10. Circumstellar dust was the obvious IR source. The lack of 10 ftm spectral features 
ruled out silicates, while the ~ 1000°K dust temperatures implied refractory material. Geisel 
et al (1976) concluded that the dust was graphite grains, though Lewis and Ney (1979) 
have argued for metallic iron grains. Clayton and Wickramasinghe (1976) have shown that 
micron-sized grains can form in the time available. 

Observations since then have shown a correlation between nova speed class and IR lumi- 
nosity: the faster the nova, the less prominent its IR excess. This was explained by Gallagher 
(1977) as a Stromgren sphere ionization process. The faster the nova (and the higher its 
bolometric luminosity) the larger its Stromgren sphere and the less able it is to produce 
dust. Bode and Evans (1982) expanded this model’s scope by considering the rise in nova 
effective temperature as the wind starts to thin. They showed that the distance from the 
star at which grains can condense increases with time. Expanding ejecta must catch up to 
the receding condensation distance before its declining density becomes too low to permit 
significant grain formation. Bode and Evans (1982) determined that novae with an initial 
rate of decline slower than 0.1 mag/day could make grains (in the context of their model) 
while faster novae could not. 
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Dinerstein (1986) found four correspondences between the positions of IRAS sources and 
those of 67 classical novae. The slow nova Sgr 1982 showed a warm dust shell, while the fast 
nova Mus 1983 did not. Dinerstein claimed that the detected infrared emission from HR Del 
(1967), FH Ser (1970) and DQ Her (1934) is inconsistent with a simple dust condensation 
plus cooling model. She suggested that strong IR lines might be responsible for most of the 
IR flux from old, metal-rich nova remnants. 

IV.S X-RAYS FROM OUTBURSTING NOVAE 

The visual brightness of a classical nova begins to decline once mass ejection subsides, 
after most ( > 70%) of the envelope has been ejected. As the envelope remnant collapses 
back onto the white dwarf, a large bolometric correction develops and most of the energy 
radiated by the white dwarf is in the EUV to X-ray band (Prialnik, Shara and Shaviv 1978, 
1979). A blackbody temperature of ~ 3 x 10 5 °K was predicted by these authors to be 
expected from classical novae for a few years after optical maximum. 

The first detections of erupting classical no»„e in the X-ray band have recently been 
reported by Ogelmann, Beuermann, and Krautter (1984) and by Ogelmann, Krautter, and 
Beuermann (1987). The three brightest erupting novae during the EXOSAT mission were 
all detected as soft (0.5 to 2.5 keV) x-ray sources. The observations of nova Vul 84 No. 1 
and nova Vul 84 No. 2 were particularly illuminating. Both were observed approximately 
four months and ten months, respectively, after outburst, and both brightened by factors of 
> 2 during this time. During this same period of time the optical light curves of both novae 
were declining significantly. In addition, nova Vul 84 No. 1 was observed by EXOSAT eight 
days after optical maximum and the nova was too faint to be detected at this time. 

Nova Muscae 1983 was clearly seen to be fading between 2-3 years after eruption. This 
short time scale conflicts with the shocked circumstellar gas model initially proposed by 
Brecher et al. (1977) in order to explain the transient soft x-ray sources. This cooling time 
is in good agreement with the constant bolometric luminosity model of a hot white dwarf 
remnant. In this model, the source remains at constant bolometric luminosity until the 
envelope hydrogen is all processed to helium. 

From the X-rays’ fluxes and energies Ogelmann, Krautter and Beuermann (1987) at- 
tempted to determine the white dwarf mass, and the mass of the nuclear processed material 
for Nova Muscae 1983 (see Figure 11). M w< i ~ 0.8-0.9A/q and ~ 10 -6 AfQ were esti- 

mated, though significant uncertainties remained due to the likely presence of an optically 
thick wind from the white dwarf remnant. A wind with M £ 10~ 7 Moyr -1 would decrease 
the photospheric effective temperatures in the energy range observed by EXOSAT and allow 
more massive white dwarfs, M w< i > 1 Mq (Ruggles and Bath, 1979; Kato, 1984). 

Effective temperatures T e and luminosities L of novae will, in future, provide stringent 
tests of the constant bolometric luminosity model. Future soft X-Ray missions with spectral 
capability will be able to measure changing T e and L as mass ejection slows and the remnant 
envelopes collapses back onto the white dwarf. 

IV.4 POSITRON ANNIHILATION GAMMA RAYS FROM NOVAE 

Positron annihilation gamma rays are potentially powerful diagnostics of novae (Clayton 
and Hoyle 1974, Leising and Clayton 1987). The latter authors have included improved 
estimates of the nucleosynthesis expected in novae, and were particularly motivated by the 
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approaching launch of the Gamma-Ray-Observatory as a potential probe of thermonuclear- 
powered hydrodynamics. It is on the observational predictions and results of this latter 
paper that I concentrate 

Leising and Clayton (1987) adopted the standard thermonuclear scenario of classical 
novae. Their numerical simulations were based on simplified models of the physics of the 
nova TNR which, nonetheless, included a sophisticated study of the generation and transfer 
of gamma ray photons out of the erupting, expanding envelope. The temperature and density 
in a nova envelope must be accurately calculated to determine the mode of annihilation 
of the positrons produced in the weak decay of proton-rich nuclei created during a TNR. 
These physical conditions determine what fraction of the positron annihilations result in 
line photons, and what fraction appear as orthopositronium continuum. A Monte Carlo 
simulation was used to study the photon propagation through the atmosphere, yielding the 
surface emissivity and energy spectrum. 

The two most important sources of annihilation line photons from positrons early in 
a nova outburst are Nitrogen- 13 and Fluorine-18. The emerging intensity of gamma rays 
depends on the transfer and escape of the emitted photons. Most of the escaping gamma 
rays emerge from the outer few scattering lengths because the atmosphere is extremely thick 
to Compton scattering. 

Unfortunately, the most crucial assumption in the Leising and Clayton (1987) simulation 
was the distribution of the positron emitters at the beginning of the evolution. The emitted 
gamma ray fluxes are extremely sensitive to whether convection can carry positron emitters 
to within a few gamma ray scattering lengths of the surface. Assuming that convection does, 
indeed, bring the positron emitters up to the surface, Leising and Clayton (1987) compared 
the calculated fluences of annihilation lines, integrated over ~4-40 minutes for various iso- 
topes, with the sensitivities of the Gamma- Ray Spectrometer on the Solar Maximum Mission 
(SMM GRS) and the Oriented Scintillation Spectrometer Experiment on the Gamma-Ray 
Observatory (GRO OSSE). Nova envelope models of 1 O~ 3 M 0 produced enough 13 TV and 
annihilation photons to be detected up to distances of several kiloparsecs during the 
first few hours of a nova outburst. 

Clearly the detection of gamma rays from outbursting novae would be an exciting con- 
firmation of existing theoretical scenarios. This is particularly true for nuclear decay lines, 
which would be powerful diagnostics of the TNR process and of the early stages of the 
outburst. Novae should be considered very high priority targets of opportunity for the 
Gamma-Ray Observatory. 


V. HIBERNATION AND ANCIENT NOVAE 
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V.l THE HIBERNATION SCENARIO OF CATACLYSMIC BINARIES 

Dozens of novae erupt in M31 every year (Arp 1956, Ciardullo et al. 1987). Because 
only one stellar system in thousands is the close white dwarf-red dwarf binary necessary to 
produce a nova, all novae must therefore recur thousands of times (Ford, 1978). It’s therefore 
of interest to ask the question: what do novae look like during the millenia between eruptions? 
The first in-depth historical studies were those done by McLaughlin (1960). McLaughlin 
concluded that there is no essential difference in the appearance of a nova before and after 
eruption. He also determined that the initial rise to maximum light must be extremely rapid, 
usually one ot two days. 
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A much more extensive study was carried out by Robinson (1975). Robinson examined 
the pre-eruption light curves of 33 novae which erupted from 1899 to 1971. Robinson’s 
results supported McLaughlin’s first conclusion, namely that there is no apparent difference 
in the brightness of a nova before and after its eruption. McLaughlin’s second conclusion 
clearly had to be modified. Most novae do show significant increases in luminosity between 
one to fifteen years before eruption. The pre- and post-nova eruption brightness of novae has 
been shown to be ~ Lq (Warner, 1987), mostly emitted by the disk and due to accretion 
energy. However, the canonical picture of old novae has recently been cast into serious doubt 
by the recoveries of two much older novae: CK Vulpeculae (nova 1670) (Shara, Moffat and 
Webbink, 1985) and WY Sagittae (nova 1783) (Shara el al. 1984). 

Once a white dwarf ejects most of its hydrogen envelope in a TNR, the remaining enve- 
lope settles back onto the white dwarf and begins to cool on a time scale of order 100 years 
(Prialnik, Shara and Shaviv, 1978). Thus CK Vul and WY Sge, as the two oldest recovered 
novae, are much more likely to be typical of “between outburst” systems than any 19th or 
20th centry novae. Both objects are amongst the intrinsically faintest known cataclysmic 
variables. WY Sagittae, at B ~ 19 has M ~ 7-8 (Shara et al. 1984, Kenyon and Berriman 
1988). CK Vul (nova 1670) has M ~ +10, i.e. a luminosity only 10 -2 Lq- Certainly in the 
latter case, mass transfer is diminished by at least two orders of magnitude relative to that 
seen in most old novae. In addition, WY Sagittae (nova 1783) and several other old novae 
have shown at least one dwarf nova outburst (Shara et al. 1984; Livio 1987). If WY Sge and 
CK Vul are typical of most old novae, then these observations suggest that 

1) old novae decline in brightness and in mass transfer rate on a time scale of a few centuries 
after eruption, and 

2) the long term evolution of cataclysmic binaries is going to be dramatically affected by this 
variable M “hibernation scenario”. 

Two other sharp disagreements between observation and theory are neatly resolved by the 
hibernation scenario. The space density of cataclysmic variables found in galactic surveys, 
of order 10 -6 pc“ 3 (summarized by Patterson, 1984) is 2-3 orders of magnitude lower than 
that deduced from simple nova theory and from the eruption frequency of novae in M31 
(Bath and Shaviv 1976). The only way cataclysmic binary space densities of 10"® pc -3 
can be compatible with the observed number density of erupting novae in the Galaxy is if 
the recurrence time of classical novae is < 10 3 years. This in turn demands high accretion 
rates, with average M > 10~ 8 AT©yr -1 onto the white dwarfs in cataclysmics. Numerical 
simulations of mass accretion onto white dwarfs demonstrate that M> 10 -1 ®M©yr _1 keeps 
the accreted envelope too hot for strong degeneracy to develop (Shaviv and Starrfield 1988, 
Prialnik et al. 1982), and, therefore, never leads to nova eruptions. Instead, continuous 
burning or mild TNRs puff the envelopes of the mass accreting white dwarfs out into red 
giants. 

A hibernation scenario resolves these apparent conflicts by hypothesizing that cata- 
clysmics hide from most surveys during most of the millenia between eruptions. They do 
this by 1) shutting off mass transfer a few centuries after eruption, 2) keeping mass trans- 
fer off for centuries or millenia, and then 3) turning mass transfer back on again for a few 
centuries before eruption. While this might seem very contrived there are sound physical 
reasons for a cataclysmic binary to behave this way. 

Before explaining why mass transfer rates should vary so dramatically in cataclysmic 
binaries, let’s step through the scenario and describe what happens at every stage of a nova 
outburst and the subsequent evolution of the underlying binary. 
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1) When a nova erupts, the underlying white dwarf’s hydrogen rich envelope suffers a TNR. 
In a very short period of time the Eddington luminosity is reached, and mass ejection 
continues for weeks or months until most of the envelope is ejected. 

2) The visual luminosity cf the star then starts to decrease as the bolometric correction 
increases. This is because the mass ejection rate from the white dwarf decreases and 
observers see progressively deeper into the atmosphere of the material still surrounding 
the white dwarf. 

3) Eventually the white dwarf stops ejecting mass, though its envelope remains hot and 
luminous for about a century after eruption. The very hot white dwarf irradiates its 
red dwarf companion, and this irradiation helps keep mass transfer high for one to two 
centuries after the nova outburst (Kovetz, Prialnik, and Shara 1988). 

4) After a few centuries the rate M of mass transfer from the red dwarf to the degenerate 
dwarf decreases. As M drops below the critical value required for dwarf nova eruptions, 
such eruptions begin. In particular, when M £ 10~® Mgyr -1 , accretion disk instabilities 
should occur and old classical novae should be able to undergo dwarf nova eruptions. At 
least eight such systems are known and Livio (1988) has recently tabulated them. It is 
not clear just how low the mass transfer rates in cataclysmic binaries become. They may 
decrease to the point that dwarf nova eruptions become infrequent or stop altogether. 

5) After a long interval of very low M (that is, hibernation) the mass transfer rate starts 
to increase; dwarf nova eruptions begin again; and if the mass transfer rate keeps on 
rising, dwarf nova eruptions will eventually be suppressed when M > 10~ 9 A/©yr -1 . 
Above this rate the system would appear as a nova-like variable, as it would have the 
spectroscopic properties of an old nova with no sign of a visible shell and no history of 
a recent TNR. 

6) Eventually a nova-like variable (that has gone through high M, then low M, then high 
again M) must erupt again as a classical nova when the pressure and temperature at 
the hydrogen-rich envelope base becomes large enough to trigger a TNR. In this way, 
a cataclysmic binary can re-cycle itself thousands of times, becoming in turn a dwarf 
nova, a hibernating system, a born-again dwarf nova, a nova-like variable, and finally a 
classical nova. 

The detailed questions we now seek to answer are: 

1. Why are post-novae bright? 

2. Why does M decrease in the centuries following eruption? 

3. How low can M in a cataclysmic binary drop? 

4. Why do the systems stay in hibernation for thousands of years before M rises again? 

5. Why are novae as bright before eruption as after? 

6. Do all cataclysmic binaries and nova systems undergo a hibernation phase? 

The answers to these 6 questions are as follows: 

1. During nova eruptions, and for decades afterward, the red dwarfs in cataclysmic binaries 
are irradiated with hundreds of times more luminosity than they themselves produce. 
The simulations of Kovetz, Prialnik and Shara (1988) followed the structure and evolu- 
tion of red dwarfs in cataclysmic binaries blasted by tremendous amounts of luminosity 
from their erupting and subsequently cooling companions. Simulations of the time de- 
pendent irradiation of three red dwarf models of 0.25, 0.5 and 0.75 Af© indicate that the 
red dwarfs swell ~ 1% in radius because of the irradiation. This swelling is maintained 
for several decades after outburst. As a result, the mass transfer rates forced by irradia- 
tion after nova eruption are found to be enhanced by up to two orders of magnitude due 
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to the irradiation. As the white dwarf cools on a time scale of a century, the enhanced 
mass tre '*r declines. 

2. M will w^rease in a cataclysmic binary if the mass-transferring red dwarf loses contact 
with its Roche lobe. This can occur in the following natural way after a nova eruption. 
During the weeks or months of nova eruption and mass ejection from the white dw-rf, 
the mass of the underlying binary system is, of course, decreasing. Of order 10 -4 -10~ 5 A/q 
is lost. This results in a weaker gravitational interaction between the two stars in the 
binary. Assuming spherically symmetric ejection from the white dwarf, Shara et al. (1986) 
investigated the effects of mass and angular momentum losses durng a mass ejection phase. 
In most cases mass loss dominates, and the separation a of the binary stars increases during 
eruption. The relationship between the separation increase Aa and amount of mass loss AM 
is given by 

Aa AM 

rsj 

a ~ M ' 

The key point is that the red and white dwarf’s separation increases as a result of the 
mass loss. The theory also predicts that the orbital period of the binary should fractionally 
increase by about 10“* to 10” 5 . Such an increase has recently been observationally confirmed 
by Schaefer and Patterson (1983). BT Mon (nova 1939) suffered a fractional period increase 
AP/P ~ +3 X 10 - ” 5 when it underwent eruption. This is exactly in the sense and size 
predicted by Shara et al (1386). 

Another possibility to detach the red dwarf from its Roche lobe is the injection via 
synchronization torques of spin angular momentum of the white dwarf into the binary orbit 
(Lamb and Melia 1987). This can occur if the timescale for synchronizing the white dwarf 
is much shorter than that for the radius of the secondary to change appreciably. 

An increase in separation Aa/a ~ 10 -4 -10 " s is much less than the increase in the size of 
the red dwarf atmosphere due to irradiation by its hot white dwarf companion. As a result, 
the red dwarf atmosphere completely fills its Roche lobe and mass transfer proceeds at an 
enhanced rate. The essential point to note is that the irradiated red dwarf swells enough 
to continue to fill its Roche lobe and to transfer mass at a high rate for a few white dwarf 
cooling times, i.e. a few centuries after the nova eruption. 

During the next phase of evolution of the binary, M decreases as the irradiation luminos- 
ity from the cooling white dwarf decreases. The red dwarf atmosphere starts to contract back 
to its pre-eruption value. Furthermore, the increased separation of the red and white dwarfs 
in the binary helps to break contact of the Roche lobe with the red dwarf atmosphere. Al- 
though the computed values and times dependence of M are very model dependent, the long 
term trend does seem to be clear: M decreases on a time scale of one or two centuries after 
eruption. Once M does decrease to the critical value lO~*A#0yr -1 , dwarf nova eruptions can 
begin. If the mass transfer rate observed in old nova systems, 10”® A/qj/t* -1 , acts for only 
a century or two after eruption then the scenario predicts that only a few xlO “®M q is ac- 
creted during this phase. During the subsequent hibernation this initial envelope will become 
rather cold and degenerate , and it will have ample time to diffuse into the carbon-oxygen or 
oxygen-neon-magnesium white dwarf. It is this rather small amount of mass, located deep 
in the C-0 or O-Ne-Mg outer layers, and accreted in the century or two after eruption that 
is crucial to triggering the next TNR millenia later. 



3. How low can M become in a cataclysmic binary? The intrinsically faintest cataclysmic 

known, CK Vul (nova 1670), must have M < (Shara, Moffat and Web- 

bink 1986) to account for its extremely low luminosity. J.ower values might lead to 
systems which appear as non-interacting, short-period white dwarf-red dwarf binaries. 

4. In an attempt to understand the behavior of accreting novae with variable M, Prialnik 
and Shara (1986) studied the effects of rapid accretion and then hibernation in the 
preoutburst history of classical novae. As in previous simulations, no hibernation time 
and M = 10 -8 M©j/r -1 produced a mild TNR, yielding a red giant rather thar * nova. 
The hibernation models produced dramatically different results. As the hibernation 

times increased, so did the amount of diffusion of hydrogen into the white dwarf carbon- 
oxygen cores. When mass transfer was turned on again at the rate of 10~ 8 Mq yr ~ 1 after 
hibernation (M = 0) times of 10 3 , 10 4 and 10 5 yr, violent TNRs occurred in each of the three 
cases. This occurred despite the fact that all accreted mass was deposited at 10 - ® Af©yr -1 . 
Much of the mass necessary to trigger the nett nova eruption arrives at the white dwarf 
surface in the ~ 10 s yr before eruption when M again becomes high, i.e., about 10 -8 M©yr -1 . 

Livio, Shankar and Truran (1988) have extended the work of Prialnik and Shara (1986) 
in studying the case of mild hibernation. They found that a reduction in the accretion 
rate by a factor of about 100 over a period louger than a few thousand years is generally 
sufficient to ensure nova-type outbursts even in the presence of high pre-outburst accretion 
rates. Thus a total shutdown of mass transfer is not necessary to ensure strong TNRs once 
mass transfer begins again. 

5. The emergence from hibernation to a dwarf nova and eventually a nova-hke state, as M 
increases, is due to the red dwarf resuming contact with iis Roche lobe. This occurs as 
gravitational radiation and angular momentum losses through magnetic stellar and disk 
winds drive the red and white dwarf together. In most cases, magnetically driven winds 
are the dominant effect and they return the stars back to their original separation on a 
time scale of ~ 10 4 years. Once M reaches the high pre-nova eruption state of about 
10 -8 Af©yr -1 , only a thousand years or so are needed until the critical envelope mass is 
achieved and the next nova eruption occurs. Warner (1987) has suggested that there is 
a natural limit to the mass accretion and mass transfer rates in cataclysmic binaries of 
about 10 -8 Af©yr -1 because so many systems have converged to this value. This might 
occur if 10~ 8 M©yr~ 1 is the upper limit to mass transfer in an accretion disk, and could 
explain why M is so similar before and after nova eruptions. 

6. Livio and Shsra (1987) studied which cataclysmic binaries should be most likely to 
undergo hibernation. These are the systems which suffer the largest reduction in mass 
transfer rat~, after a nova event, caused by the increase Aa in orbital separation. The 
assumptions used in the derivation were: 

1) The secondary fills its Roche lobe; 

2) The secondary obeys a lower main-sequence mass-radius relation; 

3) The mass ejected in the nova outburst is equal to the mass necessary to trigger a 
TNR; 

4) The Shara et al. (1986) assumptions of conservation of mass and angular momentum 
transfer during a nova outburst; 

5) The Papaloizu and Bath (1975) red dwarf scale height H (as .» function of binary 
parameters). 

These assumptions lead to a complex functional dependence of Aa/H on the orbital 
period P and the binary mass ratio q. It was found numerically that for a gi ien orbital 
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period, systems with a higher mass ratio are more likely to enter hibernation. For a given 
mass ratio the systems with shorter orbital periods are more likely to enter hibernation. 

Livio (1988) has stated that a change in the separation of the red and white dwarfs may 
not yield as large a mass transfer rate change as previously claimed. He suggested that the 
periods of high M before and after eruptions are due to irradiation of the red dwarf by the 
white dwarf. The hibernation scenario (in this viewpoint) is unchanged, but the mechanism 
to produce M is variable red dwarf irradiation rather than variable binary separation. 

The hibernation scenario explains quantitatively the apparent discrepancy between mass 
accretion rate observations and the theory of mass accretion onto pre-novae; the surprising 
rarity of nova discovered in surveys; how metallicity can vary from nova to nova; and how 
moderate and fast novae can appear even in mctal-poor globular clusters and galaxies. The 
latter effect, of course, is due to the diffusion of hydrogen into the degenerate dwarf core 
during the hibernation phase. 

The greatest advantage of the hibernation scenario is that it offers a unifying picture of 
dwarf novae, nova-like variables and classical novae. In the hibernation scenario all three are 
the same system; we simply see metamorphosis between one cataclysmic type and another 
as the mass transfer rate in the binary changes over the millenia. 

V.2 HUNTING FOR NO VAE SEEN B Y ANCIENT ASTRONOMERS 

A clear prediction of the hibernation scenario is t v at very old novae (much older than CK 
Vul) should be intrinsically much fainter than ex-novae 10-100 years after eruption. Search- 
ing for novae which erupted before the invention of the telescope is a difficult task. European 
observations of novae prior to 1600 A.D. are almost non-existent. We must therefore rely on 
records kept by Japanese, Korean and Chinese astronomers. 

For a nova to have been noted by Oriental astronomers, it must have reached at least 
third magnitude. With an amplitude of about 12 magnitudes, which is typical for most 
novae, very old novae should now be no fainter than V ~ 15 if the hibernation scenario is 
incorrect. If the hibernation scenario is correct, then these objects may be as faint as 20th 
magnitude. 

Over the last several years, I and my colleagues (M. Bode, F. R. Stephenson, M. Pot- 
ter and A. Moffat) have been conducting a search for ancient novae. We are attempting 
to recover about a dozen objects observed over the last two thousand years by Oriental 
astronomers. Deep U and B Schmidt plates of those areas of the sky that the Oriental 
astronomers reported to contain tailless, stationary “guest” stars (to exclude comets) have 
been obtained. Scanning and processing of the plates is presently underway. 

The first object we’ve searched for is the possible nova of September 19, 712 A.D. The 
astronomical treatus of the Chiu-t’ang-shu records the following: “At night the moon was 
eclipsed and it was complete, a star entered into the moon’s spirit.” The star was not a 
planet and there is no star brighter than magnitude 7.0 in the computed path of the moon 
that night (Stephenson 1986). Our U and B plates reveal no object which is brighter on 
the U plaie than the B plate down to B = 18. If the hibernation theory is wrong, and our 
ancient colleagues did indeed see a nova being occulted, we would expect to have found a 
15th magnitude nova-like variable. There clearly is not one in the 0.5 degree 2 of the sky 
where the Oriental astronomers may have recorded a nova in 712 A.D. This, of course, is 
not direct evidence that the hibernation theory is correct, but the result is consistent with 
the predictions of the hibernation scenario. If in any of the remaining survey fields we find 
bright nova-like variables, serious doubt would be cast on the hibernation scenario. If, on the 
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other hand, for a large number of fields with certain or near certain ancient novae, we can 
’ recover no nova-like variables down to 19th or 20th magnitude, support for the hibernation 

scenario of cataclysmic binaries will have been gained. 

A system which might represent a case in which hibernation did not occur is 0623+71. 
£ This object is a blue, 12th magnitude nova-like variable surrounded by a nebula (Ellis, 

| Grayson and Bond 1984). If the nebula is, in fact, a shell ejected by a nova many centuries 

, ago, then 0623+71 could represent a binary system in which accretion at a relatively high 

*1 rate persists. Alternately, the system may be a planetary nebula with a binary nucleus 

* (Bond 1988). Unfortunately there is no historical record of the outburst of this object, and 

its age remains problematical. 

i Bruhweiler and Sion (1986) have suggested that V471 Tau may have erupted as a nova 

many centuries ago. If this suggestion is correct then this detached red dwarf-white dwarf 
Hyades binary is now in a state of deep hibernation (t.e., M = 0). 

An exciting recent discovery by Bragaglia et al. (1988) is that two white dwarfs, out of 
twenty examined for radial velocity variations, have close red dwarf companions. One of the 
two (WD0034-211) is probably a weakly mass-transferring system. This could be a system 
in deep hibernation long after an eruption. If so, the implied space density of hibernating 
novae is (very roughly) .05 that of white dwarfs, or 10 -3 — 10 - * pc~ 3 — just that predicted 
by the hibernation scenario. 


VI. RECURRENT AND MAGNETIC NOVAE 

VI. 1 RECURRENT NOVAE 

The nature of the recurrent novae has been reviewed in great detail by Webbink, Livio, 
Truran and Orio (1987). This paper offers a broad overview for anyone interested in recurrent 
cataclysmic binaries. I will only touch on a few of the highlights discussed by Webbink et al. 
(1987). 

Theorists and observers generally agree that dwarf novae are accretion-powered events 
(probably due to accretion-disk instabilities quickly dumping mass onto white dwarfs), while 
classical novae are powered by thermonuclear runaways on the surfaces of white dwarfs. The 
outburst mechanism for recurrent novae is less clear. Recurrent novae suffer outbursts on 
time scales of decades to ~ one century, intermediate between the dwarf novae and classical 
novae. Erupting recurrent novae brighten by ■ — 7—1 1 magnitudes, a larger range than all 
dwarf novae but smaller than most classical novae. Finally, recurrent nova speed classes 
range from very slow (T Pyx) to extremely fast (U Sco, V394 CrA, T CrB, and RS Oph are 
amongst the fastest known novae). Webbink et al. (1987) have concluded that the “class” of 
the recurrent nova is in fact an extremely heterogeneous group of objects. This, they point 
out, is largely due to the lack of a rigorous definition as to what constitutes a recurrent nova. 
They offer the following criteria, which any object must satisfy in order to be classified as a 
recurrent nova: 

i) at least two distinct recorded outbursts, with absolute magnitudes at maximum compa- 
rable to those of classical novae ( M v £ - 5.5), and 

ii) ejection of a shell at outburst, with ejection velocity comparable to those of classical 
novae ( Veject £ 300 km s -1 ). 

These criteria cleanly distinguish recurrent novae from classical and dwarf novae, and 
also from symbiotic stars. 
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By carefully reviewing the observations of eleven objects usually referred to as recurrent 
novae, Webbink et al. were able to eliminate all but five from the class. Only two were 
considered to be caused by TNRs on the surfaces of massive white dwarfs: U Sco and T 
Pyx. Accretion-powered events probably underly the eruptions of T CrB and RS Oph. These 
latter objects are suggested to be giants transferring bursts of mass to their main-sequence 
companions. Finally, V1017 Sgr appears to be a hybrid case which exhibits both accretion 
i events and TNRs. Recently V3S4 CrA underwent a second outburst very similar to those of 

i U Sco (Sekiguchi et al. 1988). It is likely a TNR-powered nova. 

1 Webbink et al .' s (1987) main criterion for distinguishing thermonuclear from accretion 

. events in yet to be discovered recurrent novae is the spectrum at optical maximum. TNR 

, powered recurrent novae must occur on white dwarfs with masses close to the Chandresekhar 

mass, and with envelope masses that must be extremely low (10 -,i -10 -7 M©). The density of 
such envelopes at optical maximum is so low that electron scattering opacity can exceed the 
true absorption opacity by several orders of magnitude. The continuum spectrum observed 
must be formed deep within the envelope, and fc expected to depart significantly from a 
black body. Webbink et al. (1987) thus predict an emission-line spectrum, as observed in T 
Pyx, for TNR powered recurrent novae. 

In the case of accretion powered recurrent novae, the resultant shock ejection should 
produce exceptionally high velocities, decreasing with time; and the appearance of high- 
excitation coronal line emission and X-ray emission during decline. 

Finally, because of the very short accretion timescales involved in recurrent novae, T Pyx 
and Li Sco offer an outstanding opportunity to differentiate between the diffusion and the 
shear-mixing mechanisms of CNO enrichment of novae ejecta. This was discussed in §11.3. 

VI. 2 NOVAE ON STRONGLY MAGNETIC WHITE DWARFS 

The evolutionary relationship between magnetic and non-magnetic nova systems is un- 
known. The white dwarfs in the former can have strong magnetic fields. The magnetic 
cataclysmic variables divide naturally into two subclasses. The DQ Herculis stars (recently 
reviewed by Warner 1985) are inferred to have magnetic fields of 10 5 -10 7 gauss. The AM 
Herculis stars, recently reviewed by Liebert and Stockman (1985), typically have magnetic 
fields of 2-3 X 10 7 gauss. An extensive review of the evolution of magnetic cataclysmics has 
recently been published by Lamb and Melia (1987). 

DQ Herculis, the prototype of the DQ Her systems, erupted as a classical nova in 1934. 
Novae were not known to be binary systems in 1334 and we are unaware of any attempts 
to measure periods in the system then. Nevertheless, DQ Her impressed upon theorists that 
; a classical nova could indeed erupt on a white dwarf with a moderately strong magnetic 

field. Until 1987, however, it was not clear that a classical nova could erupt on a extremely 
magnetic white dwarf, i.e., in an AM Herculis binary (Livio 1983; but see Livio, Shankar 
and Truran 1988). 

Kaluzny and Semeniuk (1987) recently carried out CCD-imaging photometry of V1500 
Cygni. The fight curve of this old nova, now at 17th magnitude, revealed sharp features 
reminiscent of AM Herculis sources. Following up on this discovery Stockman, Schmidt and 
Lamb (1988) (hereafter SSL) carried out a program of polarimetric observations. Periodic 
variations of optical circular polarization with semi-amplitudes of ~ ±1.5% were found 
by SSL. The period of polarimetric variation is 1.8% shorter than the stable photometric 
period, first obtained by Patterson (1979), and refined recently by Kaluzny and Semeniuk 
(1987). SSL ascribe the circular polarization to optical cyclotron emission from the accreting 
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magnetic white dwarf primary. However, most of the system light and the photometric 
variability are due to the heated companion star. 

Nova Cygni 1975 achieved a maximum absolute magnitude M„ = —10.2 as one of the 
most luminous novae ever seen (see the recent review by Lance, McCall, and Uomoto 1988). 
A 3.3 hour variation was seen in the early nova emission lines and then in the optical 
continuum. Patterson (1978, 1979) observed the period to decrease from 0.141 days to 0.138 
days during the first year after eruption and then to increase and finally stabilize at 0.140 
days towards the end of 1977. The main theoretical thrust of SSL’s work has been to offer 
a simple, reasonable explanation for the observed period changes. 

The stable 0.140 day photometric period observed since 1977 is probably due to viewing 
the heated face of the orbiting red dwarf at moderate inclination (Patterson 1979). This 
is energetically plausible if the central white dwarf luminosity was L wc i ~ 1G 3 L© in 1978. 
Detailed models of the cooling of post-novae (Prialnik, 1986) and the resulting heating of 
the secondary (Kovetz, Prialnik and Shara 1988) support this high white dwarf luminosity, 
as well as the high resultant luminosity of the secondary a decade after eruption. A detailed 
comparison of the optical polarimetric and spectrographic observations with the theory is 
given by SSL. Here we concentrate on their explanation of the period changes due to asyn- 
chronous rotation set up in the V1500 Cyg system. The period evolution is sketched in 
Figure 13. 

The absolute visual magnitude of V1500 Cyg was > +8.3 before eruption (Wade 1987, 
Lance, McCall and tlomoto 1988). AM Her and VV Pup have M v ~ 8.6 in iheir bright 
states (Patterson 1984). These values correspond to mass transfer rates of M ~ 1O -1O M0 
per year. The angular velocity fi of the strongly magnetic white dwarf was probably phase 
locked to the orbital angular velocity fij of the binary before outburst (SSL). A massive 
white dwarf (Shara 1981, Prialnik 1986) and sufficient accretion time to permit hydrogen 
envelope enrichment in CNO isotopes are the other requirements which enabled Nova Cygni 
to become one of the fastest novae ever observed. 

SSL have proposed the following scenario for the observed period changes in V1500 Cyg. 
As the nova erupted the envelope expanded to engulf the entire binary. This significantly 
increased the moment of inertia of the white dwarf. The relative change in the moment of 
inertia is given by 


AJ _ 4x //A 4 /a\ /GMj\ 3 /? 4 / R\ 2 ( 1 \ 
T ~ 30 UJ V3/\ 4 ) 1 -0\Rg) \Mi) 


where p. is the mean molecular weight per baryon, Jb is the gas constant, a is the radiation 
constant, R is the outer radius of the envelope and R g is the radius of gyration of the 
original white dwarf. Assuming that (3 — 1 — Znuc/^£cW is — 0.1, where L n uc and 
are the nuclear and Eddington luminosities, a density p(r) ~ 10~^ gm cm -3 in the vicinity 
of the secondary, R ~ the binary separation, M\ — 1.4M© and R g ~ 7 x 10~ 4 Rq, SSL 
find A I /I ~ 10 ■ By assuming that convection, driven by the steep temperature gradient 

during the initial expansion phase dominates energy transport, then the angular velocity ft 
will also decrease by Aft/ft ~ — 10“ 2 . Thus, immediately following the TNR and resultant 
envelope expansion the white dwarf rotates more slowly than the binary. SSL identify the 
0.141 day period observed six days after outburst with the slower white dwarf rotation period. 
They also suggest that a hotspot, which illuminated the expanding nebula and which was 
seen photometrically, might have been caused by a local area of enhanced nuclear burning 
(Shara 1982), if the magnetic field or a convective plume could transport the luminosity to 
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the photosphere. The expanded envelope coupled the white dwarf to the orbiting companion 
and was quickly spun up from the 0.141 day period to the 0.140 day binary period. Some 
of the envelope was ejected, but some radiated enough thermal and potential energy to sink 
back onto the white dwarf; this caused the photometric period to decrease from 0.140 days 
to the polarimetric, white dwarf rotation period of 0.137 days. By 1977 the white dwarf 
had become faint enough in the visible part of the spectrum for the heated secondary to be 
the dominant system light source. At that point, the photometric variability stabilized at 
the orbital period of 0.140 days. SSL have suggested that it might be possible to phase the 
current polarimetric observations with the 1976 photometric data to test their model and to 
identify the dominant accretion/nuclear burning pole in V1500 Cyg. 

Clearly much more work remains to be done on this fascinating object, both by theorists 
and by observers. It is particularly gratifying to see the strong overlap and interaction 
between theory (local thermonuclear burning, strong magnetic fields, expanding envelope 
interaction with the secondary, irradiation of the secondary) and the painstaking observations 
of subtle polarization and minute changes in orbital period. Invaluable information and 
understanding is clearly to be gained by looking in detail at novae not only during the 
months after eruption but years and decades afterwards. 


VII. NOVAE IN CLUSTERS AND GALAXIES 

VII. 1 NOVAE IN GLOBULAR CLUSTERS 

Two nearly certain classical novae have been seen in globular clusters: T Sco, which 
erupted in M80 in 1860, and the nova of 1938 in M14 (Hogg and Wehlau 1964). In ad- 
dition, two dwarf novae are known in globulars: V101 in M5 (Margon, Downes and Gunn 
1981; Shara, Moffat and Potter 1987) and V4 in M30 (Margon and Downes 1983). These 
objects may be hibernating novae (see the section on hibernation). Intriguing photometric 
evidence for the presence of cataclysmics in the inn^r regions of M71 is seen in the beautiful 
color-magnitude diagram of Richer and Fahlman (1988); but spectra now indicate that the 
candidates are subdwarf O stars (Richer 1988). 

Novae which erupted in galactic globular clusters before ~1920-1930 (when systematic 
photographic surveys of globulars began), during the fall and winter months when globulars 
are less observable, and those novae which erupted near the centers of clusters with high 
central surface brightness would not have been detected. Trimble (1977) concluded that 
most globular cluster novae are in fact not observed, and that the nova rate per unit mass 
in giobular clusters is significantly higher than the rate in the Galaxy. 

Bright X-ray sources are known to be much more common in globular clusters than 
anywhere else in the galaxy (Katz 1975). This is believed to be brought about by the 
process of tidal capture in globular cluster cores (Fabian, Pringle, and Rees 1975; Sutantyo 
1975). Two unbound stars passing sufficiently close to each other in the core of a dense 
globular cluster can raise significant tides on each other and divert enough orbital energy 
into dissipative stellar oscillations to bind the pair. Succeeding periastron passages repeat 
the dissipative process, and eventually the orbit circularizes with a semi-major axis of a 
few stellar radii. This mechanism works as well, of course, for white dwarfs capturing main 
sequence stars as for neutron stars (which presumably form the bright X-ray sources). This, 
then, is a mechanism for manufacturing nova-producing binaries in globular clusters. The 
similarity of the distributions of novae and globular clusters about the central bulge of M31 
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prompted Ciardullo et al. (1986) to suggest that many or most novae are created in globular 
clusters and then ejected into the field . 

Very close binaries are believed to drive the evolution of their host globular clusters 
(Spitzer and Mathieu 1980). Because white dwarfs are much more numerous than neutron 
stars in globular clusters (and everywhere else in the Galaxy), there are probably tens or 
hundreds of times more cataclysmics in globulars than bright X-ray sources. Thus close, 
interacting binaries may control the structure and the evolution of all or most Galactic 
globular clusters. Finding these objects, and studying their properties is clearly important. 

The first optical counterpart for a globular classical nova was recently identified by 
Shara et al. (1986). Nova 1938 in M14 was identified on the basis of positional coincidence, 
brightness and blue color to be a B=20.2 object within one arc second of the outbursting 
object photographed in 1938. Spectrographic confirmation of the candidate was difficult 
because of the star’s faintness and because of the crowding of its field, but has now been 
accomplished (Shara, Moffat and Potter 1988). An even more challenging task will be the 
recovery of the nova of 1860 which erupted a few arc seconds from the core of M80. This 
object’s spectrographic confirmation will certainly require the high spatial resolution of the 
Hubble Space Telescope, and will be observed by Dr. B. Margon in the first two years 
following HST launch. 

Hertz and Grindlay (1983) iden f: fied low-luminosity Einstein X-ray sources in globulars, 
and suggested that they are short-period cataclysmic (white dwarf -f red dwarf) binaries. 
Deep CCD multicolor photometry (Margor °nd Bolte 1987) and high time resolution pho- 
tometry (Shara et al. 1988) in candidate and other fields have detected no optical cataclysmic 
binary candidates in u> Cen or in 47 Tuc. These data do not exclude Hertz and Grindlay’s 
(1983) suggestion, but raise the following question. r f the low-luminosity X-Ray sources in 
globulars are really like cataclysmic variables, then why aren’t they blue, and why don’t they 
flicker rapidly in brightness? 

Novae recur on timescales of ~ 10 4 -10® years. Even if only 2-3 novae erupt per century in 
all the globular clusters in the Galaxy, a few nova-producing binaries must still exist in most 
globular clusters in the Galaxy. Scott and Durisen OP78) pointed out that the eruptions of 
globular cluster novae would have a profound effect on the gas content of globular clusters. 

Typical globular cluster stars should lose ~ .2 A/© during evolution from the main se- 
quence to the white dwarf stage. One would expect ram pressure stripping to sweep accu- 
mulated gas out of a globular cluster during passages through the galactic plane. Between 
passages, however, ~ lOOAf© of gas should build up in a typical globular. Both 21cm and 
Ha searches of globular clusters have yielded stringent upper limits significantly below a 
hundred solar masses for neutral or ionized hydrogen (Taylor and Wood 1975; Smith, Hesser 
and Shawl 1976). The latter authors reported a conservative upper limit of 5.9Af© of ionized 
gas in NGC 6388, and a “probable” upper limit of 1.8Af©. NGC 6388 has a very large escape 
velocity (~ 75 km s -1 ) (Illingworth and Freeman 1974). Its inability to retain at least 15-20 
times more gas than is observed suggests that another depletion process exists in globulars. 

Scott and Durisen (1978) have calculated the heating effect of multiple nova eruptions 
on the interstellar gas in globular clusters. They demonstrated that even a few novae in a 
globular cluster will heat the interstellar gas to over 4 x 10® °K. Such temperatures will 
deplete globular clusters of their interstellar gas on a time scale of 10®-10® years, much less 
than the 10® year galactic plane crossing time. This is probably why globulars are so gas-free. 
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VII. 2 EXTRAGALACTIC NOVAE AND THE DISTANCE SCALE 

Novae are extremely luminous at optical wavelengths, and plentiful and'straight for ward 
to find in external galaxies. Only supernovae are much more luminous, though they are much 
rarer than novae. The fastest novae are 3 magnitudes brighter than the brightest cepheids. 

Zwicky (1936) was the first to show that there is a tight relationship between the absolute 
magnitude of a nova at maximum light and its subsequent rate of decline. McLaughlin (1945) 
and Cohen (1985) have calibrated the absolute magnitude rate of decline relation for galactic 
novae. These observations have reaffirmed that classical novae are amongst the best standard 
candles available in astronomy (both because of their brightness and because of the small 
intrinsic dispersion in the maximum magnitude rate of decline relation). 

Shara (1981) was able to show on physical grounds why the absolute magnitude- decline 
time relation works. The relevant physics used is: 

1) the steady state, constant luminosity, continuous wind ejection model of novae first pro- 
posed by Bath and Shaviv (1976) (see §IV); (this physical model assumes that radiation 
pressure drives the nova-wind); 

2) the assumption that nova visual shutoff occurs because of envelope exhaustion; 

3) the TNR hypothesis and the consequent pre-nova enveolope structure; and 

4) a simple white dwarf mass-radius relation. 

The fastest novae occur on the most massive white dwarfs because these stars require 
the lowest mass envelopes to produce a TNR, and these envelopes can be ejected the most 
quickly. Thus bright novae fade rapidly, while fainter novae fade more slowly. Only one 
assumption in the derivation of the absolute magnitude-decline time relation may limit its 
validity to any particular type of galaxy, stellar population, or epoch. 

The assumption is that metallicities in the ejected envelopes can be enriched even in 
metal poor systems. Moderate and fast novae (which demand Z > Zq) have been observed 
in the globular cluster M80 and in the Large and Small Magellanic cloud. The diffusion 
and/or shear mixing processes should, and seem to work as well on carbon-oxygen white 
dwarfs in metal poor systems as in metal rich systems. 

Another assumption is the implicit one that all novae are binaries, and that binaries 
do form in all stellar systems at all times, and that some eventually evolve to the close red 
dwarf + white dwarf pair needed for a nova. The red and white dwarfs are assumed to be 
like those in Galactic novae. This is supported by the fact that novae have been seen in 
the Galaxy and in M31 (both of type Sb), in M33 (Sc), the Large and Small Magellanic 
clouds (both irregulars) and even in M80 (a globular cluster) and in NGC 205 (dE). We 
don’t know enough about binary star formation or binary evolution to be sure that binaries 
form in galaxies at large redshifts. Despite these caveats it is important to note that novae 
in nearby galaxies obey the same absolute magnitude-decline time relationship found in the 
Milky Way, with no apparent dispersion due to metallicity or galaxy type. See Figure 14- 
Liller and Mayer (198?) have been surveying the Galaxy for erupting novae since 1983. 
The deduce a rate of 73 ± 24 yr -1 in the Milky Way, larger than that observed in M31 by 
Arp (1956). 

The pioneering observations of Arp (1956) and those of Rosino (1964, 1971) lea. t a 
nova rate of 26 ± 4 novae yr -1 in M31. The density of novae was found to increase towards 
the center of the galaxy, but there was a curious “hole” found in the center of that galaxy. 
No novae were discovered by either author within 4 arc minutes of the center of M31. Arp’s 
observations also indicate that about 1/3 of the novae in M31 are fast, 1/3 are moderate 
and 1/3 are slow speed novae. 
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Recently, Ciardullo, Ford and Jacoby (1986) and Ciardullo et al. (1987) have searched 
for novae in M31 with charge coupled device cameras (CCDs) and Ha interference filters. 
The high quantum efficiency of the CCDs allow the searches to be carried out with 1 meter 
class telescopes. More important, the dynamic range of the CCDs permits the detection of 
novae even against bright backgrounds, particularly when using Ha narrowband interference 
filters. Erupting novae are 3 to 5 magnitudes brighter in Ha than in the continuum and so 
they can be seen much longer in Ha than in the continuum (Figure 15). Ciardullo et al. 
(1987) found that the number of novae per square arc minute increases monotonically to 
the center of M31 (Figure 16), and suggested that the “hole” seen by Arp and by Rosino is 
in fact an observational artifact. The Ciardullo et al. (1987) observations indicate that the 
nova rate in M31 is more than twice that previously claimed. 

Ford and Ciardullo (1988) have shown that the time averaged B and Ha luminosity 
functions of novae are power laws. They thus have little utility as distance indicators. 
However Ha observations of galaxies remain a powerful tool for determining nova rates and 
parent populations. 

A decade-long survey of the Large and Small Magellanic clouds has been carried out by 
J. Graham (1981). Patience is required as only ~l-2 novae per year erupt in the LMC, and 
only 1 every ~ 2 years in the SMC. Graham plotted the positions of all the known novae in 
the Large and Small Magellanic clouds (see Figure 17). He found that novae are not found 
in the bars of each galaxy, implying that the low mass stars there are not yet old enough 
to have formed novae. Both the speed class distribution and the spectra of a few erupting 
Magellanic cloud novae appear to be similar to those in M31 and in the Galaxy. 

The most challenging detections of novae ever reported are those of Pritchet and van den 
Bergh (1986). During a fifteen night dark run at the Canada-France-Hawaii telescope they 
detected nine novae in the Virgo cluster. The brightest of these objects appeared at 21.75, 
and the faintest at B— 25. The light curves of most of these novae are too poorly determined 
to permit an accurate distance measure to the Virgo cluster, but a rough estimate of the 
derived Hubble constant H 0 = 69 km s -1 Mpc -1 was given. These important observations 
demonstrate that novae can be detected out to cosmologically interesting distances. They 
will surely be used by Hubble Space Telescope observers as standard candles in the Virgo 
cluster and perhaps beyond. 
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VIII. SUMMARY AND PROSPECTS FOR THE FUTURE 


An enormous amount of effort has gone into trying to understand the physical mecha- 
nisms and binary system parameters which govern nova eruptions. The accretion process, 
ignored in the early simulations of novae, is now understood to be extremely complex and 
crucial in determining the character of any outburst. Thermonuclear runaways are certainly 
not one dimensional, spherically symmetric events. Rotation, tides, magnetic fields, and 
anisotropic accretion all contribute to TNRs being local events in their early stages. If, and 
how quickly localized runaways spread over the white dwarf surface is not yet known. 

The helium, CNO, and Ne abundances of nova ejecta are powerful clues both to the 
natures of the TNRs and the previous accretion history of the white dwarf. ONeMg and even 
iron-rich white dwarfs likely exist in some nova systems. White dwarf masses approaching 
the Chandrasekhar mass seem required to explain the extraordinarily rapid recurrences of U 
Sco, V394 CrA and T Pyx. X-rays are beginning to be useful diagnostics of post-eruption 
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white dwarfs, while 7-rays could be powerful probes of the nucleosynthetic products of novae 
in the next decade. 

Dramatic changes in the orbital periods of V1500 Cygni have offered intriguing glimpses 
into the interplay between expanding envelope and underlying binary. The hibernation 
scenario of cataclysmic binaries offers a unifying picture of classical novae, dwarf novae and 
nova-like variables. These systems may undergo metamorphosis on timescales of hundreds 
to thousands of years as the mass transfer rate varies. 

Studying novae formed by the tidal capture mechanism in globular clusters might pro- 
vide valuable comparisons with field cataciysmics formed by ordinary binary evolution. A 
particularly intriguing question is: Can tidal captures form cataclysmic binaries with orbital 
periods between 2-3 hours? Either binary stellar evolution does not make these binaries, or 
they are not active as mass transfer binaries when they evolve to 2-3 hour periods. Eruptive 
binaries in globular clusters have much to teach us about cataciysmics in the held. 

Continued surveys for erupting novae in M31, the Magellanic Clouds and more distant 
galaxies will give us a better understanding of the parent populations and hence evolution- 
ary time scales of close accretion binaries. They may also contribute significantly to our 
knowledge of galaxy distances out to the Virgo cluster and beyond. As we learn more about 
novae, and understand which mechanisms and binary parameters control their outbursts, 
these powerful standard candles will be used with increasing frequency and confidence. 

Bright, naked-eye novae appear about once a decade. Much has been learned about nova 
eruptions since the spectacular appearance of Nova Cygni 1975. Astronomers will be ready 
with a large battery of new observational tools and tests when the next bright guest star 
appears in the heavens. 
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FIGURE CAPTIONS 

FIG. 1. - Heavy element (Z) profiles at the onset of the thermonuclear runaway for models 
with \f w 4 = 1.25 M©. Solid hnes correspond to L w j = 1.5 x 10 _2 L©, dashed lines to 
L w £ — 7.5 x 10 -2 Lq. The initial Z- profile is the dotted step function. The mass m is 
measured in both directions away from the initial core envelope boundary. The labels 
have the suno meanings as in Figure 2. From Kovetz and Prialnik (1985). 

FIG. 2. - The resulting accreted mass (in units of 10 s Af©) and 2 ena (%) in initial parameter 
space. B is the massive 1.25M© white dwarf, H is the hotter white dwarf (.0151©) while C 
is the cooler object (.0075L©). The suffixes 9 and 11 refer to mass accretion rates of 10 
and 10 -11 M©yr -1 , respectively. Thus the model BH9 is a 1.25M© white dwarf radiating 
.0151.0 and accreting at the rate of 10 - ® Af©yr~*. By the time its thermonuclear runaway 
starts, it has accreted 6.71 x 10~® Af© and the envelope has Z = .102. From Kovetz and 
Prialnik (1985). 
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FIG. S. • Schematic model for the white dwarf in a cataclysmic variable, illustrated for the 
case where shear turbulence is not efficient at redistributing angular momentum over the 
entire degenerate region. If shear turbulence is efficient, the core will be more nearly in 
uniform rotation, and there will be no substantial degenerate subregion in the nearly 
Keplerian envelope. The solid line represents the photosphere of the white dwarf, and 
the dashed curve, a smooth extension into an optically thin inner disk. The dot-dash 
curve illustrates the schematic form of the angular velocity fi( u>). From Durisen (1977). 

FIG. 4. - Qualitative sketches of the likely hydrodynamical behavior of the matter in a 
localised thermonuclear runaway on a white dwarf. From Shara (1982). 

FIG. ft. - Evolution of a 1.25 Mg nova model’s bolometric luminosity (solid line), nuclear 
luminosity (dashed line), and visual magnitude (dot-dashed line) through a full cycle, as 
indicated and the accretion phase preceding it. The time axis is divided into intervals 
corresponding to phases of different evolutionary time scales. From Prialnik (1986). 

FIG. 0. - History of a 1.25 Mg nova model’s envelope structure. The time axis is the same 
as in Fig. 5. The ordinate is the mu above the original WD core (m > 0) and below it 
(m < 0). The solid line labeled “M” shows the total stellar mass (M^ + M) as a function 
of time. The convective zone’s boundaries are indicated by dot-daehcd curves. Hatched 
areas within dashed boundaries denote burning shells (before the onset of convection), 
where the nuclear energy generation rate exceeds 10 s ergs g -1 s -1 . The different stages 
of mass loss and the mass of the envelope remnant are pointed out by arrows. From 
Prialnik (1986). 

FIG. 7. - Low-dispersion IUE ultraviolet spectrum of Nova CrA 1981 on 1PS1 April 12. Flux 
is in units of 10 -14 erp cm"^s -1 A -1 . From Williams et al. (1985). 

FIG. 8 . • Variation of photospheric density, pressure, radius and temperature for outflowing 
winds as a function of A My, the decline stage in visual magnitudes, for a black body 
continuum at constant luminosity. Three luminosities of 0.5 L e j, and 2 L e & are 
considered. From Bath (1978). 

FIG. 0. - Magnitude, M#, in a I? filter as a function of the wind mass loss rate, M. An 
outflow velocity of 10 8 cm/s, and luminosities of 2 L e ^ (a) and 0.5 L e{ j (b) are assumed The 
stage at which oxygen and nitrogen could be photoionized by the underlying continuum, 
and thus possibly give rise to associated emission features are listed, and indicated on 
both curves. From Bath (1978). 

FIG. 10. - The infrared development of Nova Vulpeculae 1976 (NQ Vul) is typical of that 
of moderate speed novae, showing the discontinuity in the visual light curve and the 
corresponding rise in the infrared (3.5/xm). The temperature of the dust shell is also 
included (from Bode 1981, adapted from Ney a id Hatfield 1978). 
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FIG. 11 . - The range of experimentally allowed bolometric luminosity, blackbody combina- 
tions that are compatible with the average soft X-ray counting rate (3.0 ± 0.310 -3 cts _1 ) 
of Nova Muscae 1983, 462 to 707 days after the optical maximum. The boundaries, 
indicated by solid lines, are determined by: a) limiting the luminosity to L e fa for a 1.4 
Mq white dwarf as 2 . 10 38 er <7 s _1 (top boundary); b) using the lower limits of counting 
rate, Ng, and distance (lower left boundary); c) using the upper limits of counting rate, 
Njj, and distance (upper right boundary); d) using a minimum surface radius of 5.10® 
cm for the white dwarf (lower right boundary). The central line corresponds to the mean 
estimate of counting rates, Ng , and distances. The dashed curves are evolutionary tracks 
of the central stars of planetary nebulae with core masses 0.8 and 0.9 Mq, as adopted 
from Paczynski (1971). From Ogelman, Krautter and Beuermann (1987V 

FIG. 12 . - Emerging fluxes of annihilation line photons from positrons emitted by each 
unstable nucleus, versus time, for envelope of 1 O - 4 M 0 . From Leising and Clayton (1987). 

FIG. 13. - A sketch of the orbital evolution of V1500 Cygni due to its nova eruption. From 
Stockman, Schmidt and Lamb (1988). 

FIG. 14. - The absolute magnitude M# of novae at maximum light vs. the logarithm of 
time (in days) required to decline by 3 magnitudes from maximum. From Shara (1981) 

FIG. 15. - The B and Ha light curves for a typical nova in M31. This characteristic light 
curve shows that shortly after outburst the nova is much more easily detected in Ha 
than in B. From Ford and Ciardullc (1988). 

FIG. 18a. - The positions of the 35 novae detected in the Ciardullo et al. (1987) Ha survey. 
The boxes represent their CCD fields; the ellipses reflect M31’s isophotes as measured 
by de Vaucouleurs(1958). The position of M31’s nucleus is assumed to be '>(1975) = 
0"41 m 22.2 # ; £(1975) = 41° 08' OS''. Ten of the 35 novai are within 2' of .ne nucleus. 
From Ciardullo et al. (1987). 

FIG. 10 b. - The cumulative distribution of isophotal radii for our 35 novae in M31. The solid 
line is the total amount of B light contained in the survey fields vs. the corresponding 
isophotal radius. Near the center, the novae follow the light precisely; at larger distances, 
the nova density drops off faster than the light. From Ciardullo et al. (1987). 

FIG. 17. - Distribution of Novae in the LMC. The vertical cross marks the centroid and 
its extent indicates the standard error of the mean value. The diagonal cross shows the 
center of rotation used by Feitzinger (1980). From Graham (1984). 
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